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1. INTRODUCTION 

An investigation of rarefied, internal gas dynamics with emphasis 
on shaft sealing applications was initiated on May 1, 1966 at the 
University of Tennessee in the Department o'f Mechanical and Aerospace 
Engineering. The investigation was conducted for the National Aeronautics 
and Space Administration under Research' Grant 43-001-023. This report is 
the final report of this investigation with major emphasis on the 
rarefied-ga9 viscoseal. 

The initial efforts were directed at a number of basic problems. 

The results of these investigations have previously been documented 
and only the references will be repeated here. These investigations 
concerned rarefied flow through short tubes (1,2), annuli (3,4), long 
square tubes (5), and other basic geometries (6-8). From these 
investigations a more complete understanding of rarefied flow has 
provided the basis for analysis of more complex applications. 

The major effort of this study has been directed to developing 
analytical models which can be used to predict the performance of 
viscoseals operating with a rarefied gas as the sealant. Since no 
experimental data were available which were adequate for establishing 
the validity of analytical models a significant portion of this effort 
has been devoted to the development of such data. This report details 
the analytical models which were developed and also the experimental 
procedures and results. Conclusions are presented based on the 
comparison of results from analytical models with experimental results. 

2. ANALYTICAL MODELS 

Three different analytical models have been developed in this 
investigation and will be presented in detail. The first is based on 
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Hodgson's (9) work and is called the modified Hodgson Model. The next 
model is based on the work of Boon and Tal (10) and is called the Slip 
Modified Boon and Tal Model. The third model is based on a super- 
position of the groove and annulus flow components and is called the 
Annulus-Groove Model. 

2.1 Modified Hodgson Model 

In Hodgson's (9) investigation only seals with a single thread 
start were considered. The development presented follows closely that 
of Hodgson but the analysis has been generalized to include multiple 
thread start seals and also a technique is presented for solving the 
specific molecular flow rate equation for a sealed condition. 

Hodgson chose to analyze the particular configuration where a 
smooth shaft rotates within a grooved housing. It will be shown later 
that the analysis for the grooved shaft and smooth housing is identical 
to this configuration. This being the case, Hodgson's model is not as 
restricted in this respect as it might appear. 

Hodgson considers the flow in the viscoseal to be composed of 
three basic components: (1) the pressure induced flow along the groove, 
(2) the pressure induced flow over the lands, and (3) the rotor induced 
flow in the groove. This treatment of the flow is quite common and is 
exactly the way King (11) chose to divide the flow for analytical 
considerations. Figure 1 shows a representation of the flow components, 
Qt , Qp» and Q' represent the pressure induced land flow, the rotor 
induced groove flow and the pressure induced groove flow, respectively. 
Hodgson also restricts his considerations to a seal with a single thread 


start . 
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In his development of the two pressure induced flow components, 
he takes the semi-empirical approach of Knudsen (12) in describing the 
flow throughout the entire regime from continuum to free-molecule flow. 
In his long tube work Knudsen found that his experimental data could 
be described by an equation of the form 


Q' 


[y'P + e 


, 1 + VP , 
1 + £ p J 


dP 

dV 


( 1 ) 


In the continuum limit Knudsen reduced Equation (1) to Q' = y'P 
and by equating this to the known continuum solution, he was able to 
determine Y'. In a similar manner as the free-molecule limit was 
approached he reduced Equation (1) to Q* = and by equating this 

to the known free-molecule solution he determined 8' . By considering 
the slip flow regime, Knudsen was able to determine the ratio v/£' 
by noting that Equation (1) becomes Q' = (Y'P + 8'^)-—, and equating this 
to the known continuum with slip solution. Next Knudsen determined 
the difference £'-v from a consideration of nearly free molecule flow 
and thus was able to determine both and V. He then applied experi- 
mentally determined corrections to and V so that the experimentally 
observed minimum in the Q'/AP versus P curve would be correctly predicted. 
Knudsen applied the analysis above to a long circular tube. In his 
analysis Hodgson applies the identical procedure to his treatment of 
the pressure induced flow in the groove (a long rectangular duct) and 
the pressure induced flow over the lands which he takes to be a narrow 
slit . 


2.11 Pressure Induced Flow in the Seal Groove . Hodgson assumes 
that the pressure varies continuously along the axis of the seal groove 
from a value of at the high pressure end to P^ at the other 
end. In order to determine the pressure gradient 
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gradient along the groove axis, it is necessary to relate the length 

of the groove, & , to the seal length, L. A development of a viscoseal 
8 

is shown in Figure 2. Hodgson only considered seals with a single 
thread. The development that follows is generalized to any number of 
thread starts, n . 

d 

In Figure 2 line AB is drawn perpendicular to the grooves. The 
number of turns of spiral that AB crosses is equal to the total number, 
N t , of complete turns on the seal which can be expressed as 

m _ L cos a 
t " w' 4- b' ‘ 

The length of groove per turn of spiral is TTD/cos a. The total groove 
length on the seal is 

£ = L cos a _ L^D 

T cos a * w* + b' w' + b’ * 


The length of each groove is then 

p &T Lttd 

g = N N (w' + b' ) * 
s s 

It follows that 


dP 

d& 


D[ 


dP 

LttD 

N (w' 
s 


+ b') 


N (w' + b’) 
s 

TTD 


dP 

dL * 


Hodgson indicates that the groove pressure gradient in a single threaded 


seal is 


dP = w 1 + b' dP 

dl ttd dL 
8 

from which it follows that the flow in a groove of a multi-threaded 
seal is N times the flow in the groove of a single threaded seal of 
the same groove width, land width and diameter. Since there are N g of 
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2 

these grooves, the total groove flow in the multi-threaded seal is N 

s 

time the flow in a comparable single threaded seal. The total groove 
flow is then 


1 + C.P ,p 

Qg = V( B p + c inrc7 ) dT h 


( 2 ) 


where the bracketed term is the single threaded groove flow developed 
by Hodgson. The constants B, C, C^, and C£ depend on the geometry of 
the seal and the properties of the sealant and are given in Appendix A. 
Appendix A also contains a discussion of the flow models used by 
Hodgson to obtain these constants. 

2.12 Pressure Induced Flow Over the Lands . Based on the assumption 
that the pressure varies continuously along the helical groove, Hodgson 
shows that the effective pressure gradient for the land flow is 


dP 

di 


w 


+ b' 


w 


cos 


dP 

C6 

all 


This pressure gradient applies equally well to both the single and multi- 
threaded seal. The land leakage flow for a seal of any number of threads 
is then 


1 + C 3 P dP 

«L- [D P+E TT^p) f 

4 


(3) 


where the constants D' , E, and also depend on the seal geometry 
and the properties of the sealant and are given in Appendix A along 
with a discussion of their origin. 

2.13 Rotor Induced Flow. Hodgson takes a simplified approach to 


the prediction of the rotor induced flow. In Figure 3 a groove cross 
section is shown with the rotor moving over the top of the groove. 
Hodgson develops the rotor induced flow on a molecular basis , but as he 
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polnts out a continuum approach yields the same result. It is assumed 
that the rotor induced flow is the same as the flow obtained in a long 
rectangular duct of width b' and height h in which the upper wall moves 
with velocity U cos a (the component of the circumferential velocity of 
the rotor along the axis of the groove). Rather than solve the describing 
differential equation for parallel flow, Hodgson chooses to compute 
the volume flow based on an area-weighted average velocity. This average 


velocity is 


Ub ' 


- (U cos ot)b ' + (2h + b')(0) _ 

U 2 (b ' + h) 2 (b ' + h) 


cos a. 


The total rotor induced flow is thus 


0 . = N 

q R s 2 (b ' + h) 


cos a] P = N A'P 
s 


(4) 


where is generalized to a multi-threaded seal. 

2.14 Total Seal Flow . The flow rates given by Equations (2, 3, 
and 4) are superimposed to give the total flow in the seal which can 


be expressed as 


1 + CLP 

% - - N > p + c rrcjF 1 dL 


1 + C 3 p dP 

[D ' P + E TT-cT ] dL - N S A ' P (5 > 


where flow in the direction of decreasing pressure is considered positive. 

2.15 Solution of the Modified-Hodgson Equation . Unlike Hodgson's 
basic equation. Equation (5) is applicable to a seal with any number of 
threads. There are three solutions to Equation (5) which are of 
particular interest. These three cases are: (1) the flow rate through 
a non-rotating seal, (2) the flow rate through a rotating seal, and 
(3) the pressure difference across a rotating seal when the net flow 


is zero. 
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Flow Rate in a Static Seal 
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Flow Rate in a Rotating Seal 
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Pressure Difference at Zero Net Leakage 

Since a rarefied viscoseal, in Che ideal case of a true space 
environment, will normally operate with ■ 0, it would be impossible 
to maintain a zero net flow. However, the condition of zero net flow 
is of interest as far as experimentation and comparison to continuum 
performance are concerned, and in non-space applications where =jf 0. 

The maximum pressure difference under which a seal can maintain 
a zero net flow is found by solving the modif ied-Hodgson model equation 
subject to the condition that = 0. Integrating Equation (5) subject 
to this condition one obtains 


2 fl S 

r 2P + AP,N 2 C+E 2 + C„ (2P + AP) N s C( C„ 1) 2 + C,(2P + AP) E( C, ‘ 1) 

l J s 2 2 4 4 


2P - AP 


2 + C 2 (2P - AP) 


2 + C 4 (2P - AP) 


= exp [N A’L - (N B +D’)AP]. 
s s 


(9) 


In general Equation (9) cannot be solved explicitly for AP. For the 
special case of continuum flow 2P > > AP,the continuum AP can be 
obtained from Equation (9) as 


N A'L 

AP = . (10) 

n B + D 
s 

At the other extreme, free-molecule flow. Equation (5) can be reduced 
to its free-molecule limit, 


9 HP 

Q* = - (N C + EHf - N A'P. 
N s dL a 


(10a) 


Since in the free-molecule limit P is very small, the last term in 
Equation (10a) could presumably be very small; but since efficient seal 
performance requires that the rotor induced flow be of the same order 
of magnitude as the pressure induced flow, this term is retained. 
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Solving Equation (10a) subject to the condition that ■ 0, the free- 
molecule pressure ratio becomes 

r = exp[N A'L/(N 2 C + E)]. (11) 

p s s 

The AP across the seal can be expressed as 
_ r - 1 

AP = 2P -E -— . • (12) 

r + 1 
P 

Combining Equations (11) and (12), the AP in the free-molecule regime 
is 

exp [N A'L/(N 2 C + E)] - 1 

AP = 2P 5 . (13) 

exp [N A'L/(N C + E)] + 1 
s s 

Since no explicit solution to Equation (9) can in general be 
obtained, some approximate solution technique must be employed. Since 
Newton's method of approximating roots is generally a rapidly converging 
iterative method, it is employed in solving Equation (9). Newton's 
method requires that the given relationship be differentiable. 

Equation (9) can certainly be differentiated. Another very important 
requirement is the ability to make a close initial approximation to the 
solution. This requirement is particularly important with a complicated 
relationship such as Equation (9). 

In order to make a close approximation to the roots of Equation (9), 
one needs to know as much as possible before hand about the character 
of the solution. At this point two characteristics are known: (1) the 
continuum limit. Equation (10), and (2) the free-molecule limit, 

Equation (13). The solution of Equation (9) will simply define the 
behavior over the entire range of average pressures and approach the 
continuum and free-molecule limits at the two extremes. In order to 
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make a close approximation to the value of AP which satisfies Equation 
(9) for a given P, a high average pressure, P^, is chosen initially such 
Chat the continuum solution from Equation (10) is a good approximation 
to the root of Equation (9). Using this initial approximation for AP, 
Newton's method is employed to solve Equation (9) for the AP at P^. 

An incremental decrease, AP, in the average pressure is then taken 
and the AP obtained from the previous iteration at P^ is used as the 
initial estimate of the solution for AP at P^ and the iteration process 
is repeated to obtain AP„ The entire process is repeated to obtain 
AP at P^, P^, etc. 

Essentially the above process could be continued until the entire 
spectrum of pressures had been traversed, but the process encounters 
difficulties near the intersection of the continuum and free-molecule 
asymptotes. In the near continuum region, the AP versus P solution 
of Equation (9) is fairly flat, thus making the solution at P^ ^ a 
good approximation to the root at P^. In the intersection region, 
however, the rate of change is so large that the method used above 
for the initial estimation of the solution is not sufficiently accurate. 
Two simple modifications to the above method can help to ensure a 
close approximation of the root at P . An obvious modification would 
be to reduce the step size, AP. Since all previous points P^ through 
P^ _ ^ have been determined, these points can be used to extrapolate 
to an initial approximation at P . Essentially an extrapolation of 
order n - 2 could be made, but experience has shown that a linear 
extrapolation combined with successive reductions in AP is sufficient 
to ensure convergence at P . 
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After passing through the critical region in the vicinity of the 
intersection of the asymptotes, the curve essentially assumes the 
straight line predicted by Equation (13). The solution is thus complete. 

The Fortran programs for the solution of Equation (5) for the 
three cases outlined above are presented in Reference 13. The flow rate 
in a static seal is simply a special case (Q^ = 0) of the rotating 
solution. Consequently only one program is needed for the flow rate 
solutions. 

The modif ied-Hodgson model provides a convenient means of predicting 
the performance of a given viscoseal. The limiting conditions of 
Equations (10) and (13) can conveniently be used to predict optimum 
seal geometries in the purely continuum and purely free-molecule 
regimes, respectively. In the transition regime, however, Equation (9) 
would be very difficult to use in an optimization study because of the 
time consuming solution method. 

2 . 2 Slip Modified Boon and Tal Model 

In 1959 Boon and Tal published (10) a significant analysis of the 
viscoseal in the laminar continuum regime for a constant density fluid. 
The viscoseal geometry was approximated by two flat plates, one of which 
was grooved, moving parallel to each other. The flow field was developed 
from a superposition of Couette flows along and across the grooves. 

The resulting velocity distributions were integrated to obtain the 
volumetric flow rate, and the pressure generation for zero net leakage 
conditions was developed. Subsequent investigations by Stair (14) 
using liquids as the sealant showed very good agreement with the analysis 
of Boon and Tal. The data of Hodgson and Milligan (15) using air as the 
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sealant also substantiated this analysis for laminar continuum operation 
with gases. 

The Slip Modified Boon and Tal Model has been developed by formulating 
"rarefied" corrections to the laminar continuum equations developed by 
Boon and Tal. This approach is similar to that used in previous 
investigations (3, 4, 16, 17) of internal rarefied flow that have 
indicated that a single model for gas flow through tubes and annuli, 
applicable to the flow regime extending from continuum to free molecular 
flow can be derived by the inclusion of rarefied effects on the continuum 
model. 

Consider a screw formed on a shaft located concentrically ijithin a 
cylindrical housing with a radial clearance c. The annular space is 
filled with a gas and the shaft is moving relative to the housing with 
an angular velocity, Figure 4 shows a developed view of the 
viscoseal geometry. The (x, y) axes are along and normal to the 
direction of relative motion and the (£, n) axes are parallel and 
normal to the grooves. The (x, y) and (£, n) coordinates systems are 
related by: 

£ « x cos a + y sin a 

r| = y cos a - x sin a. 

Previous investigators (10, 18) have reduced the describing partial 
differential momentum equations to the Reynolds lubrication equations. 

These analyses were based on the flat plate model of Figure 2 and 
assumed steady, isothermal, two-dimensional, laminar flow. Further, 
these analyses assumed inertia forces to be negligible in comparison to 
viscous forces and neglected body forces and end effects. The 
describing mathematical model is taken as: 
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d 2 u 

dz 2 


1 3P 

v H 


(14) 


iv 

dz 2 


1 3P 
P 3r> 


(15) 


Integration of (14) and (15) gives: 


u = b H z2 + c i z + c 2 


v = h f 2 ' + C 3 Z + V 


(16) 


(17) 


The integration constants are determined by the boundary conditions „ To 
account for the non-continuum effects, slip boundary conditions are 
introduced as derived by Kennard (19). For a monatomic gas, flowing within 
two parallel boundaries, the slip boundary condition can be written 
in the following general form: 

u I n 3 u n + GA (3^) | n 
gas z = 0 wall dz 1 z = 0 


+ 2 | 

+ 4 pr3x ; l z 


0 


+ 0(A ; , 


(18) 


where u is the velocity in the tangential x direction, z is the direction 
normal to the motion, and G is a proportionality constant. If the 
temperatures of the surfaces are assumed constant and equal to the gas 
temperature then the slip boundary condition reduces in general form to 


u = u , , + GA (- 3 ^) = 0. 

gas wall dz z 


(19) 


For the particular geometry under study, the slip boundary conditions 


are: 


Along the lands: 


du 


u I n = U cos a + GA — 3—1 „= u, 

r I z = 0 dz|z = 0 1 
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du 

u I = - GX -r-£- I , E n 

r'z=h dz 1 z = h 2 

r r 


Along the groove: 


du 


u ^ = U cos a + GX . . 

g'z=0 dz 1 z = o 

du 

u | = - GX — — ^1 = u. . 

g z = h dz z= h 4 

g g 


= u„ 


Across the lands: 


dv 

v | = - U sin a + GX -r-^l E v, 

r'z =■ 0 dz 1 z = 0 1 

dv 

v | , “ - GX-t— ^1 E V-. 

r 1 z = h dz 1 z = h 2 

r r 


Across the groove: 


dv 


v | = - U sin a + GX & I = v_ 

g z = 0 dz 1 z = 0 3 

dv 


v , = - GX — “ | = v. . 

gz=h dz 1 z = h 4 

g g 


Using the above boundary conditions the following velocity 
are determined: 

Along the lands : 

1 9P , 2 , x , A ” u l, . 

u r = 2y 3£ h r z) + ( h )z + V 


Along the grooves : 

1 9P , 2 , , . A - A . 

V5il <! - y > + ( — h — + u 3* 

8 


Across the lands: 


V ""V 

2 A 


1 3P / 2 , \ L / l. x N 

v r = 2U 3^ (Z “ h r z) + ( h )z + V 

r 


Across the grooves: 

1 3P , 2 , x . A " A 

V g = 2TI s?r - h g z) + ( — h — )z + A 

g 


components 


(20) 


( 21 ) 


( 22 ) 


( 23 ) 
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The slip velocities at the vails can now be determined, and the following 
velocity distribution equations obtained: 


1 9P „ 2 , x , ,-U cos ot s z 

u r ' 2y ^ (z " h r z + ( I + 2G\'/h r ) \T 

} + GX/h r . 1 9P 2 GX 

+ U 008 “ ( 1 + Wh r } " 2y M h r IT 

1 9P , 2 , N , ,-U cos ot „ z 

u g ■ Tv ^ (z ■ Y> + ( 1"+ ' 2 G'X/h -> r 

g g 

1 + GX/h Sp ? 

■ + “ c ° s a( i + > - 2v iz h g ir 

g g 

1 9P , 2 , \ , /U sin a s z 

v r " 2il (z ' h r z> + ( 1 + 2GX/h > ~ 

r r 


1 + GX/h 

U sin a ( x + 2 GX/h ^ “ 2y iff h r h 


r . _ 1_ 9P , 2 GX 

U ' Oil 


(24) 


(25) 


(26) 


1 3P , 2 , . . ,U sin a 

2T7 aTT (z ■ h g z) + ( ' 


) 

1 + 2GA/h ’ h_ 


, 1 + GX//h g . 1 9P 2 GX 

U sin a( l + 2GX/h ) " 2h 3ff h g h * 


(27) 


Noting that the axial velocity components are: 


u = u sin ot 

y 

v = v cos a 

y 


(28) 

(29) 


the axial flow rate components Q^, Q^g» an d Q^g "^y be determined. 
The width of the flow path for the £ land flow component is (1 - y)ttD 
and the path width for the £ groove flow is yirD. The ratio of the groove 
width to the groove plus land width is defined as y. The axial 
component of the E, coordinate land flow is : 


h h 

Qr = (1 - Y)ttD / r U dz = (1 - Y)ttD / r U sin a dz. 
* 0 y 0 


( 30 ) 


Substituting Equation (24) into Equation (30) and integrating gives 
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Q 5r = (1 - Y)ttD sin cclij (||) t h 3 <- f. - A) + U cos ). 

r 


(31) 


In a similar manner: 
Q 




D sin a[^T (||) h ! (- TT ~ T> + u cos a o 4 ! 


■2y 


g g 


g 


(32) 


Q nr - (1 - Y)ttD cos a[^) r h 3 r (-f- - £> - U sin o^} 

r 

Q ng - Yttd cos a[— (^) g h g (- — - ±) - U sin a /]. 

g 

The pressure gradients in Equations (31) through (34) may be replaced 
by the more convenient axial gradients by noting that 


(33) 

(34) 


3P 3P . 

M = 37 8in a ’ 


and 


(1 - 




,3P, 

Wg 


3P 

Tj cos 


(35) 


(36) 


From the continuity of mass, the flow across the groove in the h 
direction must equal the flow across the land in this direction. Thus 
it can be shown that 


[ 

[ 


1_ ,3Pv h 3 
2 vi Wr r 


1_ A ,3 
2m Vg g 



h 

U sin a ~ ] = 
h 

U sin a jS. ] . 


(37) 


The total flow is given by 

Q " V + \g + %T + V- (38) 

Using Equations (31) through (38), the expression for the axial flow 

is : 


Q 


1 3P_ 

12y 3y 


c 3 ttD 


[t 2 (l - y)(l +—■) + (! + 


0 1 ) (B 3 ) (Yt 2 + 


U + t 2 ) 


-e-ir 
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+ U ^ D t(l - Y) (B - 1) , 

2 l + t 2 iKi + — ) 


6 3 (l + |^) 
5517 ' l1. 


where 


8 = 

Y - 


h + c 
c 

b 

w + b 


1 + 6GX 

- (1 - Y)8 3 ( 1§^-) + Y 


1 + 


t = tan a 

The Knudsen number for the viscoseal is defined as 


( 39 ) 


K 



(40) 


where seal radial clearance, c, has been selected as the characteristic 
system dimension because of the manner of its appearance in Equation 
(39) in relation to the mean free path X. 

Equation (40) can be rewritten in terms of molecular flow rate 
rather than volumetric flow rate as 


n = - 


ap c 3 ttd 


12UKT 9y 


(1 + O 


■ 2 [t 2 (l - Y)(l + ^~) + (1 + ff^HB 3 ) 


(Yt 2 + 1 )]+ up 


R 3 n . 6GX 

t(l - Y) (8 - 1) r B (1 + Fc~ ) _ 

o l cn\ 


1 + t 


HI + ^) 


1 ]. 


(41) 


Further consideration now needs to be given to the slip coefficient 
constant of proportionality "G". Various investigators have developed 
expressions for this proportionality constant ranging from 2/3 (20) 
to (2 - f)/f (19) where f is the fraction of their tangential momentum 
which molecules give up upon striking a solid boundary. Published 
values of f by Kennard (19) range from 0.79 for air flowing over fresh 
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shellac to 1.00 for air or CO 2 over machined brass. A value of f 
equal to one is taken to be realistic for the flow of gases over 
machined surfaces. Thus the value of the slip coefficient G was 
taken to be unity. 

Knudsen number as defined in Equation (40) indicates the degree 
of rarefication since it becomes large as the mean free path, X, 
becomes large. When the Knudsen value becomes small and approaches zero, 
corresponding to continuum conditions, then Equation (30) reduces to 
the continuum, no slip, solution shown by Stair (18). 

For isothermal flow the mean free path, X, is a function of the 
pressure and gas properties. The Chapman relationship (19) can be 
used to express this dependence. 


u (2 w £) 1/2 

m 

2p 


(42) 


This relationship will be used to express the right hand side of 
Equation (41) in terms of pressure. 

2.21 Variable Mean Free Path Solution . In order to solve 
Equation (41) for the flow rate in terms of the seal length and the AP 
across the seal a numerical scheme must be employed. Before integrating 
the equation is rearranged. Let P q be the pressure at which 
X = C/6G. Then 


P - 
o 


.. kT.l/2 

6GXP 6G V(2 m * 


Now Equation (41) may be written in terms of pressure 


P^'c 3 ™ 


n = - 


X f ,[t2(1 - Y)d + ±r) + (1 + ^r)Q\yt 2 + f)J 


12ykT(l + t ) 

+ UcP o P'TTDt(l - y)(B - 1) 3 3 (1 + gp) 


2kT (1 + t ) 


[- 


gP 1 

HI + pr) 


- 1 ] 


(43) 



-19- 


where 


and 


P' + P/P 


3 (1 + W r) 

ip “ (1 - Y)p + Y 

U + p-) 


Non-dimenslonalizing the length coordinate, y, by using the seal 
length, L, i.e., y' => y/L and dy = Ldy' we can write 
n = f(P’) + g(P’) 


where 


f(P') 


2 3 

Pc UD 
o 


12ykT (1 + t )L 


[t^(l - Y)(P' + 1) 


(44) 


+ (P* + i)8 3 (Yt 2 + |)j 


g(P’) = 


cup^t(l - y)(B - l) , (P' + 1/6) 


2kT (1 + t ) 


[ 6 ' 


<K1 + 1/P') 


- P’] 


Simplifying Equation (44) by defining a new flow rate parameter n we 


obtain 


• _ T dP ' 
n = f dT* + A f g 


where 


. 6yUL 
f = c 2 P 


f = f(P',6,Y,t) - [t 2 (l - Y)(P’ + 1) + (P 1 +h 

(1 + t 2 ) B 


(45) 


6 3 (Yt 2 + i) ] 


g 3 g(P' ,6,Y,t) 


t(l - Y) (6 - 

(1 + t 2 ) 




^(1 + fr) 


- P’l 



- 20 - 


Separating variables in Equation (45) and integrating over the 
seal length gives 


1 _ f 2 f(P') 

p i " “ A f ® (p,) 


dP' 


(46) 


For a given physical condition, i.e., P^, P^jBjY.t, A^, the problem 
reduces to solving Equation (46) for n. Newton's method was used by 
letting 


p ; f(p*) 

F - 1 - / , X . AT, 

P i 


dP' = 0 


A - Ag(P') 


(47) 


Newton's method predicts successive approximations to the solution of 
Equation (45) through the following relationship: 


n + 1 


4 

n 

n 


F/F' | n 


where 

f 1 = — = / 2 -MU- dP . 

dZ PJ_ (n - A f g(P'); 

Since the most useful situation for the viscoseal as a sealing device 
would be to have the net flow rate equal to zero a similar technique 
was used to solve Equation (41) for the AP which results for n a 0. 

A widely used parameter in sealing work is defined for no net flow 
in the form of sealing coefficient, A. 


A 


r: 6yUL 
2 

APc 


(48) 


Results for the variable mean free path will be presented for both net 
leakage and for no net leakage (sealing coefficient) values . 
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2.22 Mean-Value Solution . If the mean value of seal pressure is 
used in Equation (42) to calculate a mean value of mean free path* 
then Equation (39) can be approximated by replacing the local mean 
free path by an average value. This permits the simple integration of 
Equation (39) to give 


0 1_ AP c TTD 

12,i QL (1 + t 2 ) 


r s , UcTTD 

MV + T~ 

K 2(1 + t ; 


where 


f l ( V = t a - Y)(l + 6GN k ) + (1 + 


6GN ^ 21 

yWmt 2 +±>, 


f 2 (N R ) = t(l - Y)(0 - 1)(_ 


e 3 (i + 6G \) 


»K1 + 6GN k ) 


For no net flow the sealing coefficient can be written as 


A 5 


6yUL 

2 

APc 


w 

f 2 ( V 


(50) 


Results will be presented using this mean-value type solution. 

2.23 Particle Slip Correction . In the previous section where the 
slip-modified Reynolds solution was derived, slip boundary conditions 
were applied to account for the decreased intermolecular momentum 
transport at the walls. The analysis attributed a slip velocity relative 
to the walls for all the molecules adjacent to these surfaces. For 
the case of stationary walls, this implies that every molecule, on the 
average, will possess additional flow velocity due to the slip flow 
contribution. The preponderance of evidence presented by Weber (16) , 
Kennard (19), Present (20), and Fryer (21) indicate that molecules 
whose last collision was at the wall can have no slip velocity since such 
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0 


molecules are diffusely reflected from the walls. Thus only those 
molecules coming from collisions with other molecules can possess a 
slip velocity. A correction will now be developed to account for those 
molecules which do not experience slip boundary conditions due to their 
collisions with the walls. The correction will be established by 
determining the ratio of the molecule-to-molecule collisions to the 
total number of collisions, or the sum of the molecule-to- j wall plus the 
molecule-to-molecule collisions. This ratio, a, represents the 
fraction of molecules present which experience the slip boundary condition. 

The number of molecules striking the wall per unit time and unit 
nV — 

area is (20) — where V is the mean molecular speed. If one considers 
the unwrapped viscoseal geometry of Figure 4, the seal surface area is 
approximated by 


2TTDL + 2h 


TTD 

cos a 


"w + b 


) 


where the right hand term is the surface area contributed by the groove 
side walls. Thus the number of molecule-to-wall collisions in unit 
time is 

nV __ ,, , h . 

2 cos a(w + b) 


The total number of molecule-to-molecule collisions per unit time 

nV 

and unit volume is (20) — . The volume of the viscoseal is the 

volume of the annular space plus the volume of the grooves. Thu9 the 

total number of molecule-to-molecule collisions per unit time is 

nV r _ ,2 2 N , irDL*h*b , 

— [T:L(r 2 - r x ) + (w - + b)l. 



( r 2 + r l^ ^ r 2 ” r l^ » 


Since 
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and, ; due. to the .small,., clearance, of-, the- 1 , normal- viscoseal: geometry, 


D and-r^ • r£ !3 = c; 
then tHe expression- simplifies* to 


■:3J uai' 




; -: . r! . 
r . S ft ' 


.r 


Jr/.- nV-*£-' V 31,:: ,t '• ’* *■ , 0 ’<■ n ^ 1 

Y~ 7TDL [ c + — — x ] • 


.... -IT 


w + b‘ 

•' ID") ,J »V til z. It i *- - 




Hence, the proportion of the molecules which experience slip is 

,n,-! t,l ■- r> .. L Jl*- . rG'l 1-tir. '■> € J ' ‘ -- ^ - - 1 *• • ' “ ' J 


, 0 . .=» 
. J ‘ -4-»i 


1/N k [1 + y(B - 1)] 


' fT-TiTi: - o j 77 •>', u . ,.11 j- . *.■» < -■ -u -f* 

x/n k [1 + Y(B - 1)14- {[l 1 - -,- WTT) 1 

A'» r < ►..} 3-t X" ll'-OiJv'J 81 f' * \ 


(V£; a '.., " >, 

with 

b=Jl+_£ 


N = - 
K c 


.;!•» 11 : 3 (51) 

Ct - 1 ■ 


:r. 


, ,• «) a 


and 


X = — 

O r _’U w + b 


.K / 


r?r -j 


’in 




7 - (6 - 1). 

c 


In the solution of the slip-modified Reynolds' equation, it will be 
recalled that boundary conditions were not permitted on the groove side 

> - ' j •* { v ) 1 • ( j 

walls. Thus it is more in keeping with the manner of this analysis to 

omit the molecule-to-wall collisions, which occur on}_th_ese.'. surf aces 

, /. «* w - s 

/ - 1 ‘ v Q 

and contribute the h/cos a(a + b) term in the denominator of Equation 

(51). With this restriction, the equation becomes ' ; 1 

' -- \ 

1 


a = 


,!■+ N K /2[.l + C0 , - 1)Y ] 


i ' U" 


(52) 


. ^ f j . ^ t ( 1 ’ I 1 ' » * 

■In' the' 'analysis' of slip “flow through long tubes by Weber (16), 

Fryer (21), and Milligan ‘(17) ;" the resulting flow equation could be 
manipulated into two separate terms. One of these terms is .the continuum 

/ r ’ t * 

Poiseuille flow term while the second term is the slip flow contribution 
which results from application of the slip boundary conditions. 

f 

Equivalent type terms were obtained by Milligan, Cowling, and 



Wilkerson (4) in the analysis of rarefied flow in a long concentric 
annulus. For both of these geometries, the correction for the particles 
which do not experience slip was made by multiplying .the slip flow 
contribution by 0, since it is this uncorrected term which implies that 
all the molecules possess the additional slip velocity. Thus when o is 
applied in this manner the resulting slip contribution term has been 
corrected for the particles which experience wall collisions and thus 
possess no slip velocity. When this correction is made to Equation (49) 
the following equation is obtained: 


Q 


1 AP 
12y L 


c 3 ttD 

(i + t 2 > 


W + 


UcttD 

2(1 + t 2 ) 


W* 


(53) 


where 

f 3 ( V E af l ( V + c 5 (1 ' 0) ’ 

and 

f 4 (N K ) s af 2 (N K ) + C 6 (l - O). 

= yt 2 d -y)(e 3 - D 2 + e 3 d + t 2 ) 

5 3 3 ( 1 - Y) + Y 

_ = Yt(l - Y) (3 3 ~ 1)(£ ~ 1) 

6 B 3 (l - Y) + Y 

In a similar manner to that of Equation (50 ) , the sealing 


coefficient for the mean-value slip-modified Reynolds solution with 
corrections for particles which do not experience slip is 
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Often it is more desirable to express results in terms of molecular 
flow rates rather than volumetric flow rates. Equation (49) for the 
Mean-Value Slip Modified Boon and Tal solution with no correction due 
to particles which do not experience slip is 


n _ / tt N 1/2 . c TO 

— 'o— vt) t 


'2mKT 


N r (1 + t )L 


2 W ) f l ( \ ) 


r +1 

+ ( _E _) UcTO 


r p “ 1 4KT (1 + t 4 ") 


2 v W* 


(55) 


where the relationship that 
r - 1 

AP * 2 P ( r - P - + - - 1 ) , 

P 


(56) 


has been written in terms of the pressure ratio, r , across the seal in 

P 

the right hand term. Similarly, the specific leakage flow rate with 
the correction for the particles that do not experience slip is 


n 

AP 


= - ( — H -) 1/2 (- 
v 2mKT'' v 


c 2 lTD 


N k (1 + t )L 


"> f 3 ( V 


+ 


r + 1 


r)- 


UcTTD 


x — 1 9 

p 4KT(1 + t ) 


f 4 ( V‘ 


(57) 


It should be noted that for a fixed seal geometry and a given gas 
at a specified temperature, the specific molecular flow rates of 
Equations (55) and (57) are functions only of Knudsen number, the speed 
U, and the viscoseal pressure ratio, r^. As the pressure ratio becomes 
high, then (r^ + l)/(Tp - 1) approaches unity and the specific flow 
rate is a function of the speed and the Knudsen n umb er only. 

This same technique can be used to evaluate the particle slip 


correction for the variable mean-free path solution by correcting 
Equation (43) before using a numerical scheme to solve the equation. 
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2.24 Self-Diffusion . The term diffusion, as used here, refers to 
the molecular transfer which occurs due to a concentration gradient. 

If the gas is a pure unmixed gas, then the diffusion is one of self- 
diffusion. In their work concerning low density flow of gases in a 
capillary, Pollard and Present (22) suggested that at low pressures the 
total transport can be described by the superposition of the diffusive 
transport and a drift component. In his study of flow through long 
tubes, Weber (16) applied the idea of the superposition of a diffusive 
component, a slip contribution, and the viscous component. Weber 
demonstrated that his solution has the correct limiting values for 
continuum and free molecule flows and adequately describes Knudsen's 
data for long tubes (23). Milligan experimentally verified, with 
excellent agreement, the analysis technique of Weber for rarefied flow 
in long tubes (17) . Lund and Berman (24) developed empirical relations 
for the flow and self-diffusion of gases in both long and short 
capillaries by the superposition of the diffusive and drift components. 
They developed an algebraic expression which permits the direct 
computation of the Weber diffusion coefficient at any pressure and thus 
avoided the numerical integration inherent in the Pollard and Present 
treatment of self-diffusion. Lund and Berman demonstrated the adequacy 
of their model for describing self-diffusion and flow in capillaries 
in addition to the flow between flat plates. Milligan, Cowling, and 
Wilkerson (14) extended the superposition analysis technique of Weber 
to long annuli with continuing success. 

In this section, the self diffusion flow in both the grooves and 
annular space of the viscoseal will be discussed and evaluated using 
the Pollard and Present treatment as applied by Weber, 
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Annular Space Self-Diffusion 

The self-diffusion flow in the clearance space of the rarefied- 
gas viscoseal was obtained by considering this flow to be that of a 
concentric annulus. The molecular transfer was determined by evaluating 
the net number of molecules crossing a plane normal to the annulus. 

The evaluation was done by considering separately the molecules which 
come from the outer and the inner walls with the grooved inner wall 
taken as being smooth and of diameter D as shown in Figure 1. The 
details of this derivation are presented in Reference (3) and the 
procedure used was essentially that given by Weber (16) , The following 
equations were numerically integrated to obtain the diffusion 
contribution in the annular space: 


N 


net outer wall 


, 0 R 0 tt/2 

~ 2VX fr (r 2 ~ r l> / / / (R) ‘ 


R 1 ° 


sin 0 cos 2 0(1 - e R/ \ ) dR' d0 d0 


(58) 


N . . - - 2VA ~ (r 0 - Z / / (R'). 

net inner wall dy 2 1 1 „ ; a 1 n 

R 1 0 0 

l o 


.2 (2 


R_ ir tt/2 


sin 0' cos 2 0'(1 - e“ R/,N K) dR’ d0 d0\ 


(59) 


The total self-diffusion flow in the annular space is the sum of 

Equations (58 and 59) and may be arranged into specific molecular 

flow rate form using 

dn = d_ , P__ = _1_ dP £ 1 P 

dy = dy l KT^ = KT dy KT L 


V = 


2 


,2KT 1/2 
1 iTm 


and 
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to obtain 


N 

AP 


£L_ » _ u 

An H V 


2 . 1/2 


irmKT 


(r 2 " 


r l>' 


N [SUM (O.W.) + SUM (I.W.)], 

K. 


(60) 


where SUM (O.W.) and SUM (I.W.) designates the numerical integration 
of the integrals of Equations (58 and 59), respectively. 

Groove Self-Diffusion 

The self-diffusion flow in the grooves of the shaft was obtained 
in a similar manner to that of the diffusive flow in the annular space. 
The flow in the groove was determined by considering this transport to 
be that of a long groove of rectangular cross section. This neglects 
any curvature effects of the helical groove and becomes increasingly 
in error as the groove dimensions become of the order of magnitude of 
the seal diameter. The development work of the numerical scheme is 
presented in complete detail in Reference (25). 

The groove diffusion flow involves the numerical solution of the 
following type of equation: 


VN , h b ' h °° , ... 

! // t c " 7 n) ' O 


.- e/n k 8 ) 


0 0 0 0 


+ ^ (1 - e RR/N Kg)} 5 2 dz’ (dndzdO, 
rr 


(61) 


where 



is the Knudsen number of the groove flow with the groove width, b', 
being taken as the characteristic dimension. Again 


V = 


2KT.1/2 
*■ 7Tm 


and 
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dn d,P* „ 1_ dP ^1_ AP 
d£ “ dC^KT ; " KT d£ ^KT "T * 

where i is now the groove length. The relation between the length of 
one groove and the axial seal length is 

£ a 1_ ttd L M L 

N g cos cl (w + b) sin a * 

where is the number of grooves. For a multiple grooved shaft the 
total groove diffusion flow is times the flow of a single groove. 
Thus the total specific molecular flow due to the diffusion is 


N 

£P ' “ 


N sin a 

u 

TTL 


, 1/2 


TTmKT 


[SUM (Walls)], 


(62) 


where SUM (Walls) designates the numerical integration of the integral 
terms of Equation (61) to include the total walls of the grooves. 

2.25 Composite Solution . The total specific molecular flow for 
the composite solution is obtained by adding the continuum solution 
plus the slip flow contribution after correction for the molecules 
which do not experience slip plus the self-diffusion flows in the 
annular space and the grooves. Thus for the mean-value solution 

0 

( SP> composite ’ ^ < 57 > + *«• < 60 > + E 1- < 62 > < 63 > 


It should be noted in summating the flows that the contributions of each 
equation must be evaluated at the same physical gaseous state. In 
the development of the groove self-diffusion the Knudsen number, N ^g’ 
was based on the groove width, b' , rather than the seal radial clearance, 
c. Thus the corresponding value of groove Knudsen number for a given 
clearance Knudsen number is 

\g = \ ( b r) * 


(64) 
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This same technique can be used to formulate the composite 
solution for the variable mean free path analysis and other models 
such as the annulus-groove model. 

2.3 Annulus-Groove Model . In this section a simpler model than 
the Modif ied-Hodgson Model will be developed which will lend itself to 
a less arduous solution in the transition regime and will also 
correctly predict the continuum performance. The model also predicts 
the rotor induced flow in a more rigorous manner. 

The basic model will be essentially the same as the one used by 
Hodgson and which is basically the same as the simplified screw extruder 
theory presented by Carley, et al. (26). The assumptions inherent in 
these models are: 

(1) the total flow in the seal can be treated as the super- 
position of the leakage flow in the grooves, the leakage 
flow over the lands and the rotor induced flow in the 
grooves , 

(2) the pressure varies continuously along the groove and 

is constant over the cross section of a particular groove, 

(3) the groove depth is small compared to the diameter of 
the seal, thus allowing curvature effects to be neglected 
in the groove flow development, 

and 

(4) the flow in a seal with a grooved housing is identical to 
the flow in a seal with a grooved shaft (13). 

Assumption 1 obviously neglects the convective coupling of the flow 
components. The solution without this assumption is extremely compli- 
cated for even purely continuum flow (27) . Its exclusion would 
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certainly lead to an even more complex analysis when non-continuum 
boundary conditions are applied. Since the objective is the development 
of a simplified theory, the inclusion of assumption 1 is a necessity. 

One important assumption that is usually made which is not made here is 
that the groove sidewall effects are negligible. This assumption is 
one of the prime distinctions between this analysis and the Slip- 
Modified Boon and Tal Model in which this assumption is made. 

The assumptions governing the development of the component 
flows are: 

(1) the flow is steady, constant viscosity, fully developed, 
isothermal and Newtonian with negligible body forces, 

(2) the Navier-Stokes equations with non-continuum boundary 
conditions are applicable, 

and 

(3) the non-continuum boundary conditions can be expressed as (19) 
V ^wall ^ Sn^wall 

In the development of the flow components, flow models which have been 
experimentally verified will be used and reference made to their 
verification. 

2. 31 Flow in the Groove 

Although most previous investigators have initially treated the 
groove flow in two parts, this analysis will initially treat the groove 
flow as a single flow from which the two previously mentioned groove 
flows are eventually obtained. 

Based on the assumptions stated above, the Navier-Stokes equations 
reduce to the single z-momentum equation 
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a 2 v afv = i dp 

2 .-2 y dz • 

ox ay 


(65) 


The applicable slip boundary conditions shown in Figure 5 are: 

0 (66a) 


^ Co,?' 

3x 


3v 


and 


v(b'/2,y) = - GA — (b’/2,y) 

3x 

v(x, 0) = GA — (x, 0) 

3y 


v(x, h) = - U cos a - GX-^(x, h) 

3y 


(66b) 

(66c) 

(66d) 


where advantage has been taken of the symmetry about the y-axis. The 
slip coefficient, G, is usually taken as unity as will be done here. 

It is convenient at this point to non-dimensionalize the velocity 
and the coordinates and introduce an index of rarefication. The non- 
dimensionalized variables are taken as 

x=x/(b-/2) y=7/h and u = uT os' V 

A common index of rarefication, the Knudsen number is taken as 

\k - x/h - 

Introducing these new variable into Equations (65) and (66) one 


obtains 


J2 „ -2 

3 u ,_a, 2 3 u 

2 2 ^ 2 4yU cos a dz 


(b’) 2 dP _ R' 


3x 


3y 


(67) 


where 


R' = 


(b’) 2 


4ycos a 


dP 

dz 


and the boundary conditions are 


(0, y) = 0 
u(l, y) = - 


2N 


Kh 3u 
3x 


(1, y) 


(68a) 

(68b) 
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u(x, 0) = (x, 0) (68c) 

and 

u(x, 1) = - 1 - N Kh (x, 1). (68d) 

Many techniques have been used to solve non-homogeneous problems 
of the type presented above. The least arduous is a modified form of 
the method of variation of parameters (28). This method has been 
employed successfully by Ebert and Sparrow (29) to solve Equation (65) 
with four homogeneous boundary conditions , In the case considered 
here Equation (68d) presents a non-homogeneous boundary condition , but 
the method used by Ebert and Sparrow still leads to a solution. 

The method of variation of parameters is a generalization of the 
method of separation of variables, and as such, experience gained from 
employing the latter is helpful when applying the former method. In 
the method of separation of variables, the sign of the separation 
constant is taken so that the trigonometric solution is obtained for 
the homogeneous direction. In the method of variation of parameters 
the same reasoning is used in the assumption that the solution has the 
form 

CO 

u(x, y) = E 0 n (y)(cos a^x + sin a^x) (69) 

n = 1 

where the x-direction is the homogeneous direction determined by an 
inspection of Equations (68). The function 0^ (y) is a yet unknown 
function of y. One now proceeds to determine 0 n > , and C^. 

Substituting Equation (69) into the boundary condition of Equation 


(68a) one obtains 
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Z 0 (y) [C^a cos (0) - a sin (0) ] = 0 
, n / n n 

n = 1 


which implies that = 0 and that 


u(x, y) = £ 0 (y) cos a x„ 

n = 1 

I 

Substituting Equation (70) into Equation (68b) gives 

“ | 2N vh a n 

Z 0 (y) [ cos a - sin a ] = 0 

. n n a n 

n = 1 

from which it follows that for a non-trivial solution 
2N a 

Kh n , _ 

cos a sin a = 0 

n a n 

or 


( 70 ) 


a tan a = — 
n n 2N 


Kh 


(71) 


The eigenvalues, of the eigenfunctions, cos a x , are the roots of 
the transcendental relationship of Equation (71). 

The problem now reduces to determining $ n (y). Substituting 
Equation (70) into Equation (67) yields 


E ot 2 0 (y) cos a x + (-f-) 2 Z 0"(y) 
,nn n z n 

n = 1 n = 1 


cos a x 
n 


= — — Z ft cos a x 

U i n n 

n = 1 


(72) 


where £) satisfies the Fourier series 
n 


1 = Z S7 cos a x. 

, n n 

n = 1 


( 73 ) 


The Fourier cosine series coefficient, ft , must be 

n ’ 


2 sin a 
ft = n 


sin a 


n a + sin a cos a 
n n n 


= i_ r n 

a , , 2K, sin'4 a 
n 1 + n n 


( 74 ) 


a 
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where the last step comes from Equation (71). 

In order for 0 n (y) to satisfy Equation (72), it must be that 

2a 4R'ft 

0"(y) - (-— ) 0 (y) y- • (75) 

n 3 n Ua^ 

The determination of 0 thus reduces to the solution of a non- 

n 

homogeneous second-order ordinary differential equation. The boundary 
conditions on 0 n (y) are obtained from Equations (68c) and (68d). 

From Equation (68c) , one boundary condition is 

r 

0 n (°) = \ h 0n (o) - (76a) 

Substituting Equation (73) for the unity term in Equation (68d) 

leads to the second boundary condition. 




(76b) 


The solution of Equation (75) subject to the boundary conditions 
in Equations (76) is 


-9 y -0 y R'ft 

0 (y) = (C. + C. ) e n + (C,. + C„ )e n — • 

n J lnp lnu 2 np 2nu 9 


Ua 


where 


2a 

0 = — S 
n a 


R '“„ f (1 - \ h 9 n )e " 9n - (1 + N Kh e n> 


lnp Ua 2 (1 - 0 ) 2 e ® n - (1 + ) 2 e 9n 

- Kn n - 


n 


Kh n' 


n 


(77) 


a „ (1 + N KhV 


'lnu 


(1 - N Kh e n )2e " Sn - « + 


e 


n 


R'ft 


J 2np 


= [ 


] 


n 


(1 - N Kh 9 n>- C 1 + N KhV £ 

(1 ' B Kh 0 n> 2e ' 9n - (1 + Bct n 
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8(1- N e ) 

n Kh n 


J 2nu 


(1 - N_ 0 ) 2 e~ 6n - (1 + N ,9 ) 2 e 9n 
Kh n Kh n 


By substituting Equation (77) back into Equation (70) the 
dimensionless velocity distribution becomes, after considerable 
simplification , 


o CO 

(b') dP v , 
u< *- y > ‘ 2UU CC Ta dl 1 


sin a 


n 


2N 2 

n = 1 1 + Kh sin a 


]. 


cos ot X 

[ — ] ]- 


2a y a “ 

, n _ n , , 2a y 

cosh — - tanh — sinh n 

a a 


or 


n 


1 + 


2N & 


tanh 




- 1 ] 


00 sin a 

Z - ' n 


2 [ 


n = 1 


2N 2 

1 + Kh sin a 


cos a x 
][ — — ] • 


a 


n 


n 


2a y 2N Tr , a 2a y 

. , n^ Kh n , 

sinh + cosh 


,2N_ a 2 a 4N a 2a 

ri , Kh n N 2, , , n , Kh n , n 

[1 + ) ] s i n h 1 cosh 

a a a a 


*]. 


(78) 


The volume flow rate in the groove is obtained by integrating the 
velocity distribution over the area of the groove. The volume flow 
rate thus obtained is 


V’ = 


O CO 

(b ' ) h dP y _a_ . 
2y dz n = 1 a 5 


2 

sin a 


n 


-][■ 


t anh a / a 
n 


0 J L 2N ,a 

i , Kh , 2 .. , , Kh n , , n 

n 1 + sin a 1 + tanh — 

ana a 


E] 

a a 


. 2 
sin a 


= Ub'h cosa 


n 


][- 


tanh a /a 
n 


n = 1 a 


n 1 + 


2N J 1 2N„. 

Kh , 2 .. . . Kh 

sin a 1 + 

ana 


a 


-]. (79) 


tanh 


n 


The first summation term represents the pressure induced flow in the 
groove whereas the second term represents the rotor induced flow. 


i 
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The linearity of Equation (65) has been demonstrated in the 

analysis above since the solution of Equation (79) can be shown to be 

the sum of two solutions. The first term is in agreement with the 

expression obtained by Ebert and Sparrow (29) for slip flow in a 

rectangular duct with stationary walls . Milligan and Patterson (5) 

have experimentally verified the solution of Ebert and Sparrow. The 

2 

last term is the solution of V v ■= 0 with the boundary conditions given 

in Equation (66). It is thus demonstrated that it is permissible to 

obtain the two solutions mentioned above independently and then to 

add them together to arrive at the same result as Equation (79). 

At this point it is interesting to determine when the Knudsen 

number approaches zero, the continuum flow regime. The solution for 

the eigenvalues becomes 

cos a = 0 
n 


or 


nir 

a = -rr 

n 2 


n = 1, 3, 5, 7 


With these eigenvalues and with N 

Kn 


0, V’ becomes 
g 


V' 

gc 


= 16 (b ' ) _h 
1T~V 


dP 

dz 

n 



a r ^ , nTT nil. 

5 t anh 2a - 2fl ] 
n 


8Ub'h cos a 


it 


co 

I 

n = 


— r tanh 
1 n J 


nTT 

2a 


(80) 


n 1,3, 5, 7 a a a 

2.32. Correction to Rotor Induced Flow . In the analysis above 
it was assumed that the rotor velocity, U cos a, acted at the 
top of the groove (y = h) . No attempt was made to account for 
the fact that the rotor is not located at y = h but rather at 
y = h + c. Since an exact analysis of the region above the 
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groove, h £ y £ h + c, has been shown to be very complex even for 
continuum flow (27), some sort of approximation must be made to obtain 
a simple solution. The following assumptions are made: (1) the flow 
induced by the rotor in the region h £ y <_ h + c does not contribute 
to the seal discharge and (2) an effective inducing velocity less than 
U cos a acts at y = h as a result of the clearance region h _< y <_ h + c. 
This effective inducing velocity, U e , will be determined by calculating 
the average velocity in the plane y = h of a hypothetical groove of 
depth h + c and width b with the upper boundary moving at U cos a and 
with slip boundary conditions on all surfaces. The hypothetical 
groove and associated boundary conditions, are shown in Figure 6. The 
brackets in Figure 6 indicate the portion of the boundary over which 
each boundary condition is applied. 

From Equation (78) the velocity distribution induced by the rotor is 


v r 


00 „ sin a cos a x 

/ r ^ tl . 

. , ~ [ — — 1 

n = 1 n , , Kh ,2 


1 + 


sin a 


n 


2a y 2N a 2a y 

J , n J , Kh n , n- 7 

sinh -t cosh 


[- 


2a 


2N a 2a 4N , a 

r , . / Kh n N 2, , , n Kh n 

[1 + ( ) ]sinh + cosh 

d ci ci 3 


]. 


(81) 


n 


If Equation (81) is applied to a groove of depth h + c, the resulting 
distribution is 


00 „ sin a cos a x 

• t ' n ^ r n « 

v R - u co S a E -l Jij— ]■ 

„ - 1 n + _Kh 2 „ 
a n 

2a y’ 2N irlJ a 2a y' 

, , n . Kn n , n 

sinh ; — + i — cosh T~ 

a' a a' 


2a 


2 \h' tt n 9 

[1 4- ( -j ) ] sinh cosh 


2a 4N , a 
n , Kh n 


], 
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where 


a = 


h + c 


= 


1 + c/h 

V 


Kh' h + c 1 + c/h 


y' = -X . 

y h + c 


The eigenvalues of Equation (71) are still the same because 


a_ 

2N 


?N 

Kh Kh' 


The average velocity at the plane y = h or y' = 


h + c 


is 


U = / v'(x,h/(h + c) )dx. 

6 0 


Carrying through the integration, becomes 

00 sin 2 a 

U e = 2U cos a E — -[ ——2 ] 

n - 1 CT . . Kh . 2 

n 1 + sin a 

a n 


2a 2N r , a 2a 

. , n Kh n . n 

sinh + cosh 


[r + 


2N a 
Kh n. 2 


2a 


) ] s inh 


n 


Kh n 


2a 


cosh 


or 


U = [U cos a] E 
e 1 J c 


where the is equal to twice the summation in Equation (82) 
corrected groove flow is now obtained by replacing U cos 01 in 
(79) with [U cos a]E c> The total groove flow now becomes 


V' 

g 


( b » ) 3 h dP 
2h dz 


E - U b'h cos a[E E 1 
p c u 


(82) 

(83) 
The 

Equation 

(84) 


where 
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Z = 

P 


oo 

z 

n = 1 «' 


sin 2 a 


a 


n 


tanh 


n 1 + 


2N Kh 2 

sin a 

a n 


n 


a 


a 

n 


2N__, a a 

, , Kh n , , n 

1 + tanh — 

a a 


] (85a) 


Z = 

u 


OO 

E 

n = 1 


sin 2 a 


n 


tanh 


a 


3 2N 

a J , . Kh , 2 

n 1 + sin a 

a n 


][. 


_n 

a 


2N a a 

, . Kh n , , n 

1 + tanh — 

a a 


(85b) 


and 


E = 

c 


OO 

E 

n * 1 a 


. 2 a 
sin n 


2N 


n 1 + 


Kh 


o 

sin a 


n 


2 a 2N a 
, , n , Kh n 

sinh — — + cosh 

a a 


2a 


[1 + ( -- - - - ) 2 ] sinh 


2a 


2 a 4N„. a 
n Kh n 

— r + cosh — 

a a a 


]. 
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In terms of the viscoseal geometry 


dP 

dz 


dP 

dir 

g 


N (b ' + w) 
s 

ttD 


dp 4 „ dp 

dL Sln “ dL 


( 86 ) 


and since the pressure gradient in the viscoseal is assumed to be 
constant , 


dP 

dz 


sin 



(87) 


The total groove flow is obtained by multiplying Equation (84) by N . 
Noting that N g = (TID) sin a/(b’ + w) the total groove flow becomes 

„ „ _ TrD(b') 3 hAP sin 2 a v UrrDb'h cosa sina r r /n „ x 

V g " N s V g = 2yL (b’ + w) p (b~+"V) u c* (88) 

2.3 Land Leakage Flow 

The land leakage flow is taken as flow through a long annulus. 
Milligan et al. (4) have developed an expression for slip flow in an 
annulus and have obtained excellent experimental confirmation of the 
expression. The expression derived by Milligan is 
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TTD 4 dP . 

L 128M L 


( 89 ) 


where 


a a - k'V 2 n k (1 " k ’ 2) 

A t = [(1 - K' 4 ) + 1 y -- ~r— ) + [— 


K' In K’ - N (1 - K’ 2 ) 
K 


2 n 2 


[2K f (K' 2 - 1) - 2K'(ln K’)( K ' 2 - K' + 1) - 


+ 2(1 - K ,2 )(l - K') 2 N k ], (90) 

dP 

and -jj- is the effective land pressure gradient. For a single-threaded 
L 

seal, Hodgson (9) showed based on assumption (2) at the beinning of 
this chapter that 



w’ + b' 
w’ 


cos 



(91) 


As was shown in Reference (13), Equation (91) remains valid for a 
multi-threaded seal. 


Combining Equations (89) and (91) the land leakage flow becomes 


,, _ ttD^(w' + b') cos a dP TTD tw'+ b') cos a . . 

L 128m w' L dL 128M w’L \ A 


(92) 


2. 3 Total Seal Flow and Sealing Coefficient 

The total seal flow is obtained by combining Equations (88) and 
(92). The result is 


‘ K 1 Vp + A cV 


- ~ [A I E ] 
2 u c u 
c 


( 93 ) 


7TD(b') 3 h sin 2 a 
P _ 2(w' + b') 


where 
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- 'TCP 4 (b ' + w' ) cos a 
128 w* 


and 


A = 
u 


TTDb'hc z cos a 8 in a 
6(w' + b') 


The specific molecular flow rate is 

« — i 

— = J-) 

AP AP Vr 


or 


5p yLkT [A p E p + A c A L ] " _2,. m I r E "- l 1 L "u“c~u J 


-] [a 


(94) 


c kT p 

where Equation (12) is used to eliminate AP. 

A common dimensionless viscoseal performance index is called the 
sealing coefficient and is defined by 


A - 

~ c 2 AP 

where AP is the pressure difference at zero flow. By equating V T 
zero in Equation (93), the sealing coefficient is obtained as 

A E + A A, 

A = _JL_E £ ± 

All 


(95) 

to 


(96) 


u c u 

A digital computer was used to compute the specific molecular 

flow rate from Equation (94) and the sealing coefficient from Equation 

(96). In order to evaluate the summations, £ , E , and £ , it is 

p c u 

necessary to determine the eigenvalues which satisfy Equation (71). 
Newton's method of approximating roots was employed to solve Equation 
(71) for these eigenvalues at each Knudsen number. 
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3. EXPERIMENTAL INVESTIGATION 

Since no reliable performance data were available for viscoseals 
having a rarefied gas as the sealant it was necessary to conduct an 
experimental investigation. The object was to obtain performance data 

which could be used- t-o- evaluate the theoretical models. These - ~ 

investigations were conducted on multiple groove two-inch diameter 
viscoseals operating over a wide range of shaft speeds and gas densities. 
Data were obtained with no leakage through the seal which permitted the 
evaluation of sealing coefficients. In addition data were obtained for 
a wide range of net leakage conditions. All rarefied data were obtained 
using argon as the sealant. Continuum sealing coefficient data were 
obtained using air as the sealant. 

3.1 Viscoseal Test Section 

The experimental apparatus was designed to investigate viscoseal 
performance in the gas flow regime between continuum and free molecule flow. 

The viscoseal test section. Figure 7, consists of an outer housing 
with its associated vacuum pumping system surrounding a rotating grooved 
shaft. The shaft is a hollow eight-inch cantilever extension of a high 
speed spindle shaft which is belt-driven through an intermediate spindle 
by a direct current motor. The drive system is capable of seal shaft 
speeds from zero to 35,000 rpm. The speed control for the motor is 
self-regulating and maintains a selected speed within +0.1 percent. 

Since the entire test section operates under vacuum, a rubbing contact 
graphite ring seal is provided where the rotating shaft penetrates the 
housing. The ring seal is a series B-103032, Type E, manufactured by 
the Cleveland Graphite Bronze Division of the Clevite Corporation. 

Cooling water and air are supplied at the ring seal end of the spindle 
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to minimize thermal growth due to the rubbing seal friction under 
dynamic conditions. A Conflat vacuum flange is provided between the 
housing and its support for purposes of sealing and two provide minor 
adjustment capabilities for housing-shaft alignment. Other vacuum 
seals are obtained by the use of "0" rings in addition to vacuum 
sealants for threaded connections. The maintenance of a high quality 
vacuum system with essentially no atmospheric leakage, except at the 
shaft seal, was assured by the frequent use of a helium leak detector 
throughout the experimental program. 

Experimental data were obtained for two different viscoseal 
geometries consisting of a grooved shaft inside a smooth housing. 
Pertinent specifications of the rarefied viscoseals are contained in 
T ab le I . 


Table I. Dimensional Specifications of Viscoseals (All Dimensions are in 
Inches) 


Parameter 

Seal No. 

1 


Seal No. 2 


Housing diameter (+ 0.0003) 

2.0088 



2.0088 


Shaft diameter (;f 0.0001) 

2.0005 



2.0005 


"Cold" radial clearance 
(+ 0.0002) 

0.00418 



0.00418 


Axial length (+ 0.005) 

4.530 



4.530 


Groove axial width 

0.03111 

+ 

0.0003 

0.1691 + 

0.012 

Land axial width 

0.03235 

+ 

0.0003 

0.2016 + 

0.012 

Groove depth 

0.03065 

+ 

0.0003 

0.0144 + 

0.0005 

Groove helix angle 

9.30° 



19.474° 


Aspect ratio, (b cos a) /h 

1.002 



11.071 




-45- 


A schematic diagram of the overall experimental test apparatus 
is shown in Figure 8. 

3. 2 Vacuum Pumping System 

It was necessary to develop a vacuum pumping system capable of 
providing the desired low pressures and gas flow rates. This task was 
accomplished by connecting two independent vacuum pumping systems in 
parallel to a common reservoir, Figure 8. One of these pumping systems 
consisted of a single stage rotary oil-sealed mechanical pump mated with 
a three-stage high vacuum oil diffusion pump. This part of the pumping 
system was connected to the reservoir with a large diameter open-or-close 
vacuum valve. The second pumping system was composed of another single 
stage rotary oil-sealed mechanical pump mated with a positive displacement 
roots blower type vacuum pump. This portion of the pumping system was 
connected to the reservoir througha throttlable high vacuum valve. 

Through manipulation of the two connecting valves it was possible to 
regulate the pumping speed over a suitable range of downstream test section 
pressures. The complete vacuum pumping system was capable of attaining 
pressures to 10 millimeters of mercury with a blanked-off system. 

3 . 3 Instrumentation 

The radial alignment of the shaft within the housing was determined 
by the use of five proximity detectors manufactured by the Bentley- 
Nevada Corporation. These probes were located near either end of 
;>e viscoseal section. When viewed from the shaft drive end of the 
test section, the probes of the first set were located at 12, 3, 
and 6 o'clock positions, and the probes of the second set were 
located at 12 and 3 o'clock positions. The arrangement of the 
probes permitted the shaft to be aligned within the housing to a 
value of eccentricity ratio (shaft centerline deviation/mean radial 
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clearance) on the order of 0.05. The probes were also used to determine 
the thermal growth of the shaft caused by the graphite ring seal 

friction when the shaft was rotating. In addition, the probes 

i 

permitted the measurement of all vibrational movements of the shaft. 

The output gains of the detector amplifiers were individually adjusted 
and calibrated to ensure a linear output voltage of the probes as a 
function of the clearance gap. 

Pressure measurements upstream and downstream of the viscoseal 
were obtained using both McLeod gauges and absolute aneroid type 
gauges. The McLeod gauge is normally considered to be a primary standard 
(30), but after previous attempts (31) to use thermocouple gauges, 
ionization gauges and cold cathode gauges, it was apparent that McLeod 
gauges were the only instruments capable of giving the desired accuracy 
and reproducibility. The McLeod gauges utilized were the GM-100A 
gauge manufactured by the Consolidated Vacuum Corporation and are 
described in detail in Reference (31). 

Temperature measurements of the argon entering the high pressure 
end of the viscoseal test section were made early in the experimental 
program using a thermocouple. A comparison of the measured gas 
temperature with the ambient (room) temperature revealed that the error 
involved in using room temperature instead of the actual gas temperature 
was less than one percent. Therefore the gas temperature was taken to 
be 297°K in all data calculations. 

Speed measurements of the viscoseal shaft were made using a 
magnetic pickup located near the attachment nut for the drive belt 
pulley on the high speed spindle. The pulses from the pickup were 
registered using an electronic counter to obtain shaft rpm„ The speed 
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measurements were independently verified using a calibrated hand held 

tachometer in the low rpm range and a calibrated strobotac over the 

entire speed range of 0 to 35,000 rpm. 

The molecular flow rate through the viscoseal test section was 

determined using a variation of a constant pressure method developed 

by J. R. Downing (30) for the measurement of pumping speeds of vacuum 

pumps. The technique consists of applying the perfect gas equation 

of state to known gas volumes at two different time periods. In brief, 

the flow measurement is obtained by adjusting the indexed valve of 

Figure 8 to a desired setting and waiting for steady state to be 

achieved in the test section as indication by pressure and proximity 

probe measurements. The shut-off valve leading from the constant 

pressure gas reservoir is then closed and the rate of rise of manometer 

fluid in the right hand side of the manometer is observed. With 

knowledge of the manometer cross sectional area, the rate of fluid 

rise, the value of the calibrated volume, the initial and final pressures 

in the calibrated volume, and the gas temperature, the molecular flow 

rate is calculated in the following manner. 

Consider the flow measurement diagram of Figure 9. Let 

V q = calibrated volume including tank and manometer down to 

the zero deflection line, "0" 

P q = atmospheric pressure 

Hj = initial manometer deflection 

H£ = final manometer deflection 

A = cross section area of manometer tube 
m 

W = specific weight ofmanometer fluid 
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Frora Quid statics: 


P • = P + 2WH 1 

1 o 1 

P = P - 2WH„ . 

2 o 2 

The corresponding gas volumes are: 

V = V + HA 

1 o 1 m 

V = V - H A „ 

2 o 2 m 

Applying the perfect gas equation of state, 
PV = NKT, 


(97) 


on a molecular basis to determine the number of molecules within the 
system at time (1) and time (2), then. 



n 2 = 


P V 
11 

KT, 


P V 
2 2 

KT„ 


(98) 


Assuming isothermal conditions and substituting for the pressures and 
volumes. Equation (98) simplifies to, 


AN = 


H 1 + H 2 
KT 


{ A [P^ + 2W(H. - HQ] + 2WV } . 
mo 1 Z o 1 


(99) 


The flow measurement system was operated experimentally such that 



Letting 

H = H ! + H 2 

and At be the time period for the total deflection H, then the molecular 
flow rate is 


N = 


AN 

At 


H 


. (A P + 2WV } . 

AtKT mo o 


( 100 ) 


The basic premise for the measurement system rests in the fact 
that while the pressure in the calibrated volume does change slightly, 
less than + 1%, it still remains very close to atmospheric pressure. 
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Thus , the change in flow rate through the valve is negligibly affected 
by the slight decrease in the pressure within the calibrated volume 
during the course of a measurement. 

3.4 Experimental Procedure 

Experimental data of two different types were obtained for the 
operation of the viscoseals in the rarefied regime. These may be 
described as net flow leakage data and sealing coefficient data. The 
operational procedure utilized in gathering the data differed only in 
respect to shutting off the argon gas supply when obtaining sealing 
coefficient data. 

All data were acquired by the coordinated adjustment of the gas 
supply rate, the vacuum pumping speed, and the seal shaft drive speed 
followed by a time period sufficient to obtain a steady-state condition. 
Following achievement of steady state, final readings of the gas flow 
rate, pressures upstream and downstream of the viscoseal, and the 
proximity probes were made. The proximity probes were used to record 
the seal clearance which decreases slightly with increasing shaft 
speed due to the friction of the graphite ring vacuum seal- 

3.5 Data Treatment 

The viscoseal performance data for both the net flow leakage tests 
and the sealing coefficient tests were correlated versus an index of the 
flow rarefi cation, Knudsen number. The mean radial clearance, c, was 
selected as the viscoseal characteristic dimension- Thus using the 
01 paan relationship for the mean free path, the experimental values 
for Knudsen number were determined using the following equation. 


\ - = - 


;r t ) 1/2 


c 


2P c 


( 101 ) 
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The mean free path was based on the average seal section pressure 


P 


P T + P b 

2 


( 102 ) 


The molecular flow rates for the net flow leakage data were 
calculated using Equation (100). These data were revised to specific 
molecular flow rates by dividing by the pressure drop, AP, across the 
viscoseal, or 


' H [A P + 2WV ] 
N __ L m o o 

AP “ (P„ - P )AtKT 
i b 


(103) 


The sealing coefficient data were calculated from the zero net 


leakage tests using 


A , , 6PUL _ (104) 

APc (P T - P b )c 

I 

3.6 Experimental Uncertainty . A detailed analysis was performed 

to estimate the propagated uncertainty in the final performance parameters 

based on the individual uncertainties of all measured variables. The 

analysis was performed using standard statistical techniques after 

postulating a normal distribution in the uncertainty of each variable. 

This method of analysis assumes all errors are random in nature. 

Considerable effort was expended during the experimental investigation 

to minimize any error of a systematic nature. 

The propagated uncertainties, in general, are not constant and vary 

with the degree of rarefication. The uncertainty in the specific 

molecular flow rate, N/AP, of Equation (103) is + 5.1 percent at an 

Inverse Knudsen Number of 27 where the calibrated volume, V , 

o 

contribution is 84 percent of the total. The uncertainty in the 
specific molecular flow rate increases to + 8.9 percent at an Inverse 
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Knudsen Number of 0.50 where the uncertainty in the differential 
pressure, AP, contributes 72 percent of the total. The ability to 
determine the Inverse Knudsen Number of Equation (101) is estimated 

at + 4.8 percent for 1/N^ equal to 27 and increases to + 5.7 percent for 
1/N equal to 0.50. The propagated uncertainty estimate in the sealing 
coefficient, A, of Equation (104) is insensitive to the degree of 
rarefication and remains at essentially +9.6 percent. The uncertainty 
in the radial clerance, c, contributes approximately 99 percent of 
the propagated uncertainty estimate for A. It should be noted that 
the clearance uncertainty is controlled by the ability to establish 
the diameters of the seal shaft and housing and not by any limitations 
of the proximity probe system. 

All uncertainties are stated to a confidence level of 95 percent. 

3. 7 Investigations Conducted 

Sealing coefficient performance data and net leakage performance 
data were obtained for each of the two geometries. Sealing coefficient 
data were obtained over a wide range of sealant gas rarefication for 
shaft speeds of zero, 5000, 10,000, and 30,000 rpm. Net leakage data 
were also obtained over comparable ranges of rarefication and shaft 
speeds. These data are presented and discussed in Section 4 and the 
reduced data are tabulated in Appendix B. 

4. RESULTS 

* Two types of results will be presented for each of the three 
analytical models. Net leakage and sealing coefficients will be given 
as a function of the degree of rarefication. The inverse Knudsen 
number is used as the parameter indicating the degree of rarefication. 
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4 . 1 Net Leakage 

The net leakage flow occurs from the high pressure end to the 
low pressure end of the viscoseal when the viscous pumping 
action is insufficient to overcome the pressure induced flow in the 
grooves and through the annular clearance space. The net leakage data 
are shown in Figures 10 through 15. These data are presented in the 
form of specific leakage rate (units of molecules per sec-yHg) versus 
Inverse Knudsen Number. The experimental data for viscoseal Seal No. 1 
are shown in Figures 10, 11, and 12. The theoretical specific leakage 
values as predicted by the Modified Hodgson Model are shown in 
Figure 10. These theoretical values were obtained by using Equation (6) 
for the zero rpm case and Equation (8) for the 5000; 10,000; and 30,000 
rpm cases. The Slip Modified Boon and Tal Model was used to determine 
the theoretical values shown on Figure 11. Two solutions for this model 
are presented. The mean-value values were calcualted using Equation 
(55). The variable-mean-free-path values were obtained by solving 
Equation (46). In Figure 12 theoretical values for the Annulus-Groove 
Model are presented. Equation (94) was used to determine the theoretical 
values shown in Figure 12. 

The experimental data for Viscoseal No. 2 are presented in Figures 
13, 14, and 15. These figures contain data for one additional speed, 

2500 rpm. The theoretical specific leakage values shown on these 
figures were determined in the same manner as described above for 

* 

Viscoseal No. 1. 

In Section 2.25 a composite type solution was discussed for the 
Slip Modified Boon and Tal Model. This type of composite solution is 



obtained by the superposition of the corrected continuum-slip flow 
contribution with the self-diffusion flows in the annular space and 
the grooves. The corrected continuum-slip flow is determined by 
taking into account the fact that all the molecules cannot experience 

t he slip boundary condition. This procedure is ex plaine d in detail 

in Section 2.23. A detailed procedure for calculating the self-diffusion 
flows is given in Section 2.24. Tables II through VII present composite 
solution results and the uncorrected continuum-slip results. In each 
table the first column on the left is Inverse Knudsen Number. Moving 
to the right in the tables the second column is the uncorrected 
continuum— slip solution. The third column is the particle— slip 
corrected solution. The fourth column is the self— diffusion through the 
annular clearance space , the fifth column is the self— diffusion through 
the grooves, and the last column is the composite solution. 

Table II shows specific leakage results from the Slip-Modified 
Boon and Tal Model for Seal No. 1 at zero rpm. Table III is the same 
except that these results are for 5000 rpm. Tables IV and V contain 
specific leakage as predicted by the Boon and Tal Model for Seal No. 2 
at zero and 2500 rpm, respectively. Tables VI and VII contain results 
for the Annulus— Groove Model. Table VI is for Seal No. 2 at zero rpm 
and Table VII is for the same seal at 2500 rpm. These are typical 
composite solution results. The particular results shown in Tables II 
through VII were chosen because they were the cases which gave the. 
largest deviations between the composite solution and the simpler non- 
corrected continuum- slip solution. 
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Net Leakage Solutions x 10 for Seal No. 1 at Zero rpm 
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Net Leakage Solutions x 10 for Seal No. 1 at 5000 rpm 
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Table IV, Slip Modified Boon and Tal Model 
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Net Leakage Solutions x 10 for Seal No. 2 at Zero rpm 
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Net Leakage Solutions x 10 for Seal No. 2 at 2500 rpm 


- 57 - 

























o 





















rH 





















Pm 

cn 




















0) 










vO 

cn 











+ 

O 

vO 

m 

m 

O 

05 


m 

O 

p. 

vO 


00 

o 

rH 

05 

rH 



•rH 


05 

i — 1 

CM 

CM 

00 

vO 

<d 


o 

r- 

m 

cn 

rH 

m 


CM 

CM 



0) 

CM 

• 

• 

• 

• 

• 


• 

• 

• 

« 

• 

• 

• 

• 

• 

• 

0 



O 


r-H 

CM 

CM 

CM 

» — I 

1 — 1 

rH 

o 

o 

O 

rH 

CM 

cn 

n 

rH 

cn 

p 

CO 


& 

+ 

1 

1 

i 

1 

1 

1 

1 

i 

+ 






CM 

m 

p 

VO 


o 

r-H 




















o 



















c 





















o 





















-H 





















4-» 










































O 


g 

















OO 

cn 

cr 


o 











CM 

vO 

vO 

00 

cn 

1 — 1 

-d 

cn 

w 


*rH 







p. 

cn 

vO 

1 — 1 

P-. 

vO 

m 

-d 

cn 

rH 

O 

o 

Q> 

cn 


rH 

p- 

OV 

cn 

o 

-d 

cn 

rH 

1 — 1 

O 

o 

o 

o 

o 

o 

o 

o 

1 

> 

2 



<r 

CM 

CM 

1 — 1 

o 

o 

O 

o 

O 

o 

o 

o 

o 

O 

o 

o 


O 

M-J 

CO 

• 

• 

• 

» 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

g 

O 

MH 


o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

u 

•rH 



















*H 

o 

a 



















4-> 





















G 





















< 1 





















o 


c 



















CO 


o 

















m 

o 


<n 

tH 













CM 

o 

p 

o 

oo 

VO 

(U 

6 

0) 


m 

CM 

00 

vO 

05 

f — 1 

rH 

cn 

p* 

cn 

rH 

05 

00 

m 

CM 

o 

O 

•u 

i — I 

6 



m 

<M 

H 


m 

<r 

CM 

rH 

rH 

rH 

o 

o 

o 

o 

o 

O 

*H 

3 

M-l 

CM 

• 

• 

• 

• 

• 

• 

• 

• 

• 

« 

• 


• 

• 

• 

• 

• 

cn 

C 

M-l 


rH 

rH 

tH 

rH 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

q 

-H 



















a. 


n 



















u 





















O 





















o 

Q. 





















*H 

• 



















P-. 

* — 1 

u 



















<J 

CO 

u 










m 

vO 










o 


VO 

m 

CM 


05 

m 

»H 

o 

p- 

cn 





05 

00 

o 


+ 

o 

r-H 

on 

i-H 

00 

vO 

vO 

CM 

05 

o 

i — i 

vO 


CM 

o 


cn 

CM 

CM 


• 

• 


-d 

~d 

cn 

cn 

CM 

CM 

i — i 

rH 

o 

O 

rH 

CM 

cn 

in 

rH 

cn 

P 


4J 

-U 


1 

I 

l 

i 

| 

1 

\ 

1 

i 

+ 





CM 

m 

P 


c 





















O 

Cfl 




















u 

P. 





















u 





















u 





















o 









05 

vO 

vO 










CJ 


vO 

m 

CM 

o 

P-. 

o 

00 

CM 

m 

CM 

rH 

O 

05 

00 

o 

oo 

o 


p* 



CM 

CM 

CM 

CM 

o 

oo 

m 

oo 

o 

p 

in 

cn 

o 


-d 

CM 

CM 


•rH 

• 


• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


rH 

4J 


CM 

CM 

CM 

CM 

CM 

1 — 1 

*H 

O 

o 

O 

rH 

CM 

cn 

m 

rH 

cn 

P 

a, 

CO 

>-l 


| 

1 

1 

1 

| 

1 

| 

| 

i 

+ 





CM 

m 

P 

<i 

+ 

cU 





















p. 




















-U 





















c 

• 




















o 

o 


CM 

CM 

CM 

t-H 

t — 1 

rH 













o 

z 


1 

1 

| 

| 

1 

| 

O 

o 

o 

o 

O 

O 

o 

rH 

rH 

rH 

CM 





O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

O 

o 

O 





> — 1 

r — 1 

t — 1 

H 

r-H 

rH 

> — 1 

rH 

rH 

rH 

i — 1 

rH 

rH 

rH 

rH 

rH 

rH 























z 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 



1 — 1 



CO 

p* 

r-H 

cn 

r- 

*H 

CM 

cn 


m 

vO 

P 

* — 1 

cn 

P 

rH 



Table VI. Annulus -Groove Model 
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4 . 2 Sealing Coefficient 

The sealing coefficient as defined in Equation (95) is a common 
dimensionless viscoseal performance index. The relationship which 
exists between the net leakage flow and the sealing coefficient can 
best be shown by reference to Figure 16. A family of net leakage 
curves as a function of the seal pressure ratio for a fixed shaft speed 
are shown. The curves progress to the left for increasing values of 
pressure ratio until the (r + l)/(rp - 1) factor approaches unity. 

When the family of net leakage flow curves are extended to the zero 
net flow condition, a corresponding sealing coefficient condition may 
be indicated for each of the pressure ratio curves. Operating the 
rarefied viscoseal at "negative" net leakage is simply another way of 
saying that the viscoseal has now become a positive flow pump. 

In Figures 17 through 22 sealing coefficient data and theoretical 
values are presented. Figures 17, 18, and 19 show results for 
Viscoseal No. 1 at 5000; 10,000; and 30,000 rpm. Figure 17 presents a 
comparison of experimental data and theoretical values as predicted 
by the Modified Hodgson Model. Figure 18 shows the same type of 
comparison for the Slip Modified Boon and Tal Model. Theoretical 
values are shown for both the mean-value and variable mean free path 
solutions. Figure 19 is a similar comparison for the Annulus-Groove 
Model. 

Figures 20 through 22 show experimental and theoretical values for 
Viscoseal No. 2. Again data and analytical values are shown for each 
of the three models at speeds of 5000; 10,000; and 30,000 rpm. 
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5. CONCLUSIONS 

Examination of Figures 10, 11, and 12, for Seal No. 1, indicates 
that the Annulus -Groove Model gives the best agreement with the 
experimental data. The Slip Modified Boon and Tal Model gives the 
poorest agreement. All three models predict the general characteristic 
for this low aspect (groove width to groove depth ratio) ratio, 1.002, 
seal. It is concluded that none of the theoretical models are completely 
adequate for predicting the leakage through low aspect ratio viscoseals. 
The Annulus -Groove Model is recommended for low aspect ratio viscoseals 
and Inverse Knudsen Numbers greater than 5. 

For Seal No. 2, the results shown in Figures 13, 14, and 15 
indicate that the Slip Modified Boon and Tal Model gives a better 
correlation with experimental data than either of the other two models. 
The only exception is for static (zero rpm) conditions having Knudsen 
Numbers greater than 0.5 where the Modified Hodgson Model shows better 
agreement. The aspect ratio of this geometry was 11.071 as compared to 
1.002 for Seal No. 1. As anticipated the Slip Modified Boon and Tal 
Model is much better for the higher aspect geometry than for low aspect 
ratio geometries. The variable mean free path solution is only slightly 
better than the mean-value solution to the Slip Modified Boon and Tal 
Model. In view of the added complexity of the variable mean free path 
solution, it is recommended that the mean-value Slip Modified Boon 
and Tal Model be used to predict viscoseal leakage for high 
aspect ratio geometries provided the Inverse Knudsen Numbers are 
greater than 5. It is of interest to note that the 2500 rpm data for 
Viscoseal No. 2 seem to approach a constant specific leakage value with 
increasing Knudsen number. 
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In Tables II through VII a comparison between the continuum-slip 
and composite type solutions are shown for the Slip Modified Boon and 
Tal Model and the Annulus-Groove Model. These particular cases were 
chosen for presentation since they are the ones which gave maximum 
deviation between the two solution techniques. Examination of Table II 
shows that the diffusion flow contributions in both the annular space 
and the sixteen grooves are only one percent of the total composite 
flow at an Inverse Knudsen Number of 7.0 and decrease rapidly with 
decreased rarefication. At an Inverse Knudsen Number of 0.01 it can be 
observed that the diffusion contributes 92 percent of the total composite 
solution and thus the transport is essentially free molecular. This 
same trend is evident for the other two static cases as shown in 
Tables IV and VI. In Tables III, V, and VII the negative values of 
specific leakage simply indicate that the viscoseal is acting as a 
pump, i„e., the viscous pumping flow due to the grooves is greater than 
leakage flow. 

For the Slip Modified Boon and Tal Model the maximum deviation 
between the composite solution and the continuum solution is less than 
three percent for Inverse Knudsen Numbers greater than 5.0. The deviation 
increases with rarefication and has a maximum value of approximately 
twenty-five percent. This maximum occurs at 1/N = 0.01 and zero rpm 

for Seal No. 2. 

For the Annulus-Groove Model the maximum deviation is less than 

two percent for Inverse Knudsen Numbers greater than five. The maximum 

deviation was approximately sixty percent at 1/N = 0.01 and 2500 rpm 

K. 


for Seal No, 2. 
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The good agreement between the rather complex composite solutions 
and the continuum-slip solutions for Inverse Knudsen Number greater than 
five is apparent. It is concluded that the complex composite solutions 
cannot be justified for Inverse Knudsen Numbers greater than 5.0. 

Consideration of the results shown in Figures 17 through 22 
reveals several important aspects of the three analytical models. Only 
the Modified Hodgson Model and the variable mean free path solution to 
the Slip Modified Boon and Tal Model predict the general sealing 
characteristics which are evidenced by the experimental data. In 
general the models predict better sealing performance (smaller sealing 
coefficients) than what actually existed. It may be observed that the 
experimental data indicate that the sealing coefficient parameter 
becomes speed sensitive as the sealant gas becomes rarefied. This 
sensitivity to shaft speed for the rarefied-gas viscoseal is unique 
from the operation in the laminar continuum regime where the sealing 
coefficient is independent of shaft speed (14, 15). This continuum 
characteristic may be seen by noting the asymptotic behavior of the 
experimental data for 5000 and 10,000 rpm as the continuum (1/N > 

100) regime is approached. These data may be compared since the mean 
radial clearance is essentially the same for both speeds. 

Examination of Figures 18 and 21 reveals that the mean-value 
solution to the Slip Modified Boon and Tal Model can predict a seal 
pressure differential, AP, which may exceed the magnitude of the upstream 
pressure. This solution would imply a negative absolute value for the 
downstream pressure, a condition which cannot occur in reality. A 
limiting criteria to be applied to this theoretical solution to avoid 
this condition is predictable by considering the variables Involved. 
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As stated previously, the Knudsen number 


_ X P<2*/T) 1/2 


\ = 7 


2Pc 


is based on the average pressure in the seal 


P T + P 2P - AP 
P = _2 L_ = — I 


As the downstream pressure approaches zero, then AP approaches the 
upstream pressure and for this limiting condition 

— Ap 

P 0/ —L 

a, 2 * 

Thus, the Inverse Knudsen Number variable may be written as 


(1/ V 


APc 


y(2TT^ T) 1/2 


The Sealing Coefficient variable is defined as 


A 

APc 

Solving for Ap in the Inverse Knudsen Number relationship and 
substituting this into the above equation gives the limiting functional 
relationship 


A limit 


6UL 

c(27r !o T) 1/2 


'(l/\) 


(105) 


For a given seal geometry, gas, temperature, and operating - speed , the 
limiting seal coefficient function varies inversely with the Inverse 
Knudsen Number, 1/N . On a logarithmic plot such as Figure 17, this 
relation will plot as a straight line with a slope of minus unity. 
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Points which fall above and to the right of this line are possible while 
those which fall below and to the left are not possible in that a 
negative absolute pressure downstream is implied. The limiting 
relation of Equation (105) is speed sensitive and thus the limit 
condition can be established for each selected shaft speed. Limit curves 
for speeds of 5000, 10,000, and 30,000 rpm from Equation (105) are 
shown on Figures 18, 19, 21, and 22. The experimental sealing 
coefficient data agree well with these limit curves and all points 
fall within the "possible" regions. In addition, the speed sensitivity 
of the experimental data agrees with that indicated by Equation (105). 

It should be noted that the limiting function of Equation (105) 
results from a limit condition of the variables involved and applies 
to any and all viscoseal analyses using these variables. 

The limiting relationship just established is shown on Figures 18 
and 21 to indicate that the mean-value solution to the Slip Modified 
Boon and Tal Model when used in conjunction with Equation (105) is a 
useful technique for predicting sealing coefficient performance. This 
limiting relationship can also be used with the Annulus-Groove Model. 

The devaiation of predicted performance from experimental results 
is approximately the same for each of the three analytical models 
when the limiting relationship is used on the mean-value Slip Modified 
Boon and Tal Model and the Annulus-Groove Model. Theoretical values 
are not presented for the composite type solutions since there is 
negligible difference between the composite solutions for sealing 
coefficient and the continuum-slip solutions. 

Since most practical seal geometries will result in a high 
aspect ratio, near 10, it is concluded that the simplicity of the mean- 



- 66 - 


value solution to the Slip Modified Boon and Tal Model in conjunction 

i 

with the limiting relationship justifies its use. For low aspect ratio 
geometries the Annulus-Groove Model in conjunction with the limiting 
relationship is recommended. The complexity of the variable mean free 
path solution (Slip Modified Boon and Tal Model) cannot be justified 
in view of the small improvement in predicted value. The composite 
solutions are not practical for the same reasons. 

6. RECOMMENDATIONS 

In general it is recommended that the mean-value solution to the 
Slip Modified Boon and Tal Model, Equation (55), be used to predict 
viscoseal leakage for rarefied conditions resulting in Inverse Knudsen 
Numbers greater than 5. No satisfactory analytical model has been 
developed in this study to predict leakage for Inverse Knudsen Numbers 
less than 5. 

To predict sealing performance it is recommended, in general, 
that the mean-value solution to the Slip Modified Boon and Tal Model, 
Equation (50), in conjunction with the limiting relationship, 

Equation (105), be used. 

It should be noted that the variable mean-free-path solution to 
the Slip Modified Boon and Tal Model is more accurate for all conditions. 
If use is made of a high 9peed digital computer its complexity becomes 
less important and thus its use may be justified. 
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Figure 1. Basic Elements of the Viscoseal. 
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APPENDIX A. CONSTANTS FOR HODGSON'S MODEL 
The following expressions for the coefficients in Hodgson's 
Model are taken from reference (9). 


A = 


B = 


N7Tdb 2 h 


120 (b + h) 
bh 3 (b + w)K, 
12yid 


cos a 


8K l( bh) (b 4- w) 

3TTd<;b + h) C 2TTm ; 

n* = TTdc 3 (b + w) 

12Pw 


8K 2 c irdCb + w) /2 

E _ (__) 


C 1 " 


K 2 (b + h) 2 h 


kT 1 /? 

8K l bM( t^ ) [2K l b “ K 3 (b + h) 3 


Co = 


2K 1 b 


2 K 3 (b + h) 1 


C 3 = 


8 >*2 (2K 2 - «<i5 )1/2 


C 4 - 2K 2 C 3 


K = 3 ( 1 . in (a + (1 + a 2 ) 1/2 )+ a 2 In — ■ ( . 1 - a . ) 

1 8a 2 a 


2 . 1/2 


, 1 , 3 . 2,3/2.. 

+ -j (1 + a - (1 + a ) )] 


K„ = 1 - 0.63 r tanh 
3 b 2h 
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Values of K 2 
w/c 0.1 

K 2 0.036 

When w/c > 10, K 2 


0.2 0.4 0.8 

0.068 0.13 0.22 



12 3 

0.26 0.40 0.52 


4 5 10 

0.60 0.67 0.94 


The groove flow coefficients, B, C, Cp and C 2 are obtained 
from a consideration of flow in a long rectangular tube. The continuum 
groove flow coefficient, B, is obtained from the continuum solution for 
Poiseuille flow in a long rectangular tube. From a consideration of 
purely molecular flow, the f ree-molecule flow coefficient, C, can be 
derived. The basic f ree-molecule flow equation used by Hodgson is 
given in Reference (32) and is written as 


16 AT, I/ 2 

3 n'2iTnr S dl * 


(106) 


Equation (106) is attributed to Knudsen (32, p. 35), but the constant 
K , which is equal to K~ , is obtained by Clausing (32, p. 40). 0 and 

S are the groove cross sectional area and perimeter, respectively. 

The two remaining groove flow constants, and C 2 , are obtained in 
the same way Knudsen determined v and £ in Equation (1) . As was 
pointed out in Section 2.1, the determination of the ratio v/£ requires 
knowledge of the continuum with slip solution. Hodgson does not take 
an exact approach at this point since he determines C-^/C,^ from the 
continuum with slip solution for flow between infinite parallel flat 
plates rather than flow in a long rectangular tube. 

Hodgson obtains the land leakage flow coefficients, D, E, C^, 
and Cp from a consideration of flow in a thin slit-like tube. In this 
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case the continuum coefficient, D, is obtained from the continuum 
solution for Poiseuille flow between parallel flat plates. The same 
basic model, Equation (106), for f ree-molecule flow is used to determine 
the free-molecule flow coefficient, E. now becomes , which is 

also attributed to Clausing. 0 and S now represent the cross sectional 
area of the land leakage passageway and its perimeter, respectively. 

The constants and are also determined using the same approach 
Knudsen used to determine V and £ in Equation (1) . For the continuum 
with slip solution, Hodgson once again uses the continuum with slip 
solution for flow between infinite parallel flat plates which in this 
case is the proper solution to use. 
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Net Leakage Reduced Data for Seal No. 1, Zero rpm - Argon Gas 


Point No. 

Mhg 

P b 

yhg 

1/n k 

N/AP 
x 10' 

1 

- 16,800 

7230 

24.68 

18.16 

2 

19,150 

7705 

27.57 

20.30 

3 

21,400 

8155 

30.32 

21.60 

4 

22,600 

8430 

31.87 

22.08 

5 

17,000 

7465 

25.09 

17.45 

6 

12,150 

565 

13.05 

10.60 

7 

9,370 

538 

10.15 

9.06 

8 

8,100 

522 

8.84 

8.24 

9 

6,850 

512 

7.54 

7.28 

10 

5,450 

500 

6.10 

6.34 

11 

4,180 

495 

4.79 

5.50 

12 

4,200 

490 

4.80 

5.46 

13 

2,950 

487 

3.52 

4.65 

14 

1,870 

482 

2.41 

3.81 

15 

13,800 

1250 

15.40 

11.31 

16 

1,740 

255 

2.10 

3.47 

17 

1,700 

260 

2.01 

3.60 

18 

1,720 

275 

2.05 

3.59 

19 

752 

227 

1.00 

3.23 

20 

428 

132 

0.57 

3.02 

21 

318 

129 

0.46 

2.83 

22 

251 

131 

0.39 

2.40 

23 

216 

131 

0.36 

2.48 

24 

4,000 

131 

4.23 

5.26 

25 

2,320 

172 

2.55 

4.15 

26 

1,460 

166 

1.66 

3.46 

27 

1,025 

168 

1,22 

3.30 

28 

600 

172 

0.77 

2.77 

29 

475 

131 

0.62 

3.20 
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Net Leakage Reduced Data for Seal No. 1, Dynamic Speeds - Argon Gas 








Data 

Point 

No. 

P T 

yh g 

p b 

Phg 

1/N 

K 

N/AP 
x 10“ U 

c 

in. x 10 3 


10,000 rpm, 


c x 10 3 


3.87 in. , r 
1_R 


146.4 


1 

9,380 

335 

8.83 

0.558 

3.81 

2 

6,410 

60 

6.03 

0.117 

3.80 

3 

13,200 

86 

12.9 

1.50 

3.90 

4 

16,700 

93 

15.9 

2.74 

3.87 

5 

28,700 

150 

27.5 

7.80 

3.87 

6 

32,500 

166 

30.9 

9.80 

3.87 

7 

8,050 

64 

7.78 

0.237 

3.91 

10 

31,000 

233 

30. 7 

10.15 

4.01 


5,000 

rpm c x 

10 3 = 3.85 in. , 

r = 176.6 
P 


11 

16 , 800 

63 

16.6 

7.29 

4.01 

12 

11,200 

42 

10.6 

4.14 

3.84 

13 

7,900 

33 

7.18 

2.61 

3.69 

14 

6,150 

30 

5.88 

1.88 

3.88 

15 

3,700 

27 

3.51 

0.953 

3.84 

16 

1,950 

26 

2.03 

0.540 

3.76 

17 

1,250 

29 

1.32 

0.290 

3.94 


30,000 

rpm, c x 

10 3 = 3.29 in. 

, r =807 
P 


18 

44,000 

53 

36.0 

0.451 

3.33 

19 

49,900 

64 

40.9 

1.36 

3.33 

20 

39,800 

49 

31.4 

0.038 

3.20 



Rarefied Sealing Coefficient Reduced Data for Seal No. 1 - Argon Gas 

■ '.t-i Tii' g"i •m.-miiincsO 



Data Point 
No. 

P T 

V h 8 

P b 

yhg 

1/H K 

A 

Re"- c ' -T J7 

c . 

° in. x 10 


o> 



10. 

,000 rpm, 

c x 10^ = 

3.98 in. 



SP1 

7,160 

267 

7.47 

40.2 

1.01 

4.18 

SP2 

9,860 

1220 

11.1 

31.8 

1.49 

4.15 

SP3 

5,200 

275 

5.01 

68.3 

0.67 

3.80 

SP9 

14,200 

3495 

16.9 

28.3 

2.28 

3.98 

SP10 

17,200 

6605 

22.6 

29.3 

3,02 

3.94 

SP13 

19,200 

7755 

25.5 

27.1 

3.42 

3.94 

SP14 

8,700 

44 

8.29 

36.4 

1.10 

3.91 

SP19 

7,350 

44 

6.96 

43.1 

0.93 

3.91 


5000 rpm, c 

x 10 ^ » 3„ 

99 in. 



SP4 

2,600 

166 

2.69 

60.8 

0.18 

4.04 

SP5 

11,000 

4730 

15.6 

22.5 

1.03 

4.11 

SP6 

3,880 

176 

3.94 

39.9 

0.26 

4.04 

SP7 

3,160 

40 7 

3,47 

53.7 

0.23 

4.04 

SP15 

2 ,850 

50 

2.72 

56.7 

0.18 

3.90 

SP16 

13,300 

8245 

20.6 

30.2 

1.39 

3.98 

SP17 

3,920 

146 

3.69 

46.5 

0.25 

3.73 

SP18 

4,400 

382 

4.65 

36.6 

0.31 

4.04 

SP20* 

4,200 

559 

4.47 

43.6 

0.30 

3.90 

SP21* 

7,000 

1835 

8.44 

29.6 

0.57 

3.97 

SP22* 

7,550 

1710 

8.94 

25.5 

0.60 

4.01 

SP23* 

13,100 

8145 

20.7 

29,8 

1.39 

4.04 

*Argon blanket 

on graphite ring seal. 




30,000 rpm, c x 10^ 

=> 3.41 in. 



SP11 

35,200 

104 

29. 3 

34.7 

11.8 

3.44 

SP12 

38,100 

68 

30.9 

33.5 

12.5 

3.37 
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Continuum Sealing Coefficient Tabulation Viscoseal No. 1, Gas-Air 



Data 

Point 

No. 

Speed 

rpm 

AP 

in.-H 2 0 

A 

R 

e 

c 

c 3 

in, x lCr 


Test Date 6-25-69 


2 

4,905 

2.45 

23.7 

88.1 

3 

8,080 

3.74 

25.6 

145.4 

4 

11,500 

5.15 

26.5 

207.3 

5 

13,200 

5.74 

27.2 

238.1 

6 

19,600 

8.49 

27.4 

354.7 

7 

23,400 

10.38 

26.6 

424.4 

8 

25,400 

11.45 

26.2 

461,3 

9 

29,600 

14.27 

24.5 

539.5 

10 

31,800 

16.50 

22.8 

581.2 

11 

34,600 

19.58 

20.9 

634.9 

12 

21,600 

9.45 

27.0 

391.4 

13 

9,700 

4.45 

25.8 

174.6 

15 

4,800 

2.40 

23.7 

87.1 

16 

6,400 

3.19 

24.5 

116.2 

17 

7,900 

3.69 

25.4 

143.5 

18 

9,700 

4.40 

26.0 

176.4 

19 

11,600 

5.12 

26.8 

211.1 

20 

13,150 

5.74 

27.1 

239.5 



Test 

Date 7-2-69 


22 

1,272 

0.73 

20.6 

23.1 

23 

3,219 

1.69 

22.6 

58.5 

24 

4,644 

2. 34 

23.5 

• 84.5 

25 

6,399 

3.08 

25.0 

116.5 

26 

7.869 

3.68 

25.2 

143.4 

27 

9,582 

4. 35 

26.0 

174.7 

28 

11,562 

5.13 

26.7 

211.0 

29 

13,269 

5.83 

27.0 

242.5 

30 

7,392 

3.51 

24.9 

134.6 

31 

1,278 

0.74 

20.4 

23.2 





- 100 - 


Contlnuum Sealing Coefficient Tabulation Viscoseal No. 1, Gas-Air 



Data 

Point 

No. 

Speed 

rpm 

AP 

in.-H 2 0 

A 

R 

e 

c 

C 3 

in. X 10 J 


Teat Date 9-2-69 


32 

1,566 

0.87 

21.5 

28.1 

33 

3,540 

1.80 

23.3 

63.7 

34 

5,130 

2.57 

23.6 

92.4 

35 

6,690 

3.17 

25.0 

120.6 

36 

7,980 

3.70 

25.5 

144.0 

37 

9,750 

4.40 

26.3 

176.0 

38 

9,990 

4.50 

26.3 

180.4 

39 

13,800 

6.01 

27.2 

249.7 

41 

14,100 

6.21 

26.9 

255.2 

42 

19,230 

8.66 

26.3 

349.1 

43 

21,180 

9.34 

26.8 

384.8 

44 

23,400 

10.41 

26.6 

425.7 

45 

5,040 

2.50 

23.9 

90.8 

46 

1,278 

0.71 

21.3 

23.0 

47 

3,540 

1.81 

23.2 

63.7 

48 

5,010 

2.46 

24.1 

90.2 

49 

6,570 

3.12 

25.0 

118.4 

50 

7,740 

3.60 

25.5 

139.6 

51 

9,720 

4.35 

26.4 

175.5 

52 

11,490 

5.04 

27.0 

207.6 

53 

13,830 

5.99 

27.4 

250.2 

54 

2,580 

1.34 

22.8 

46.4 

55 

510 

0.28 

21.6 

9.2 



-101 


Continuum Sealing Coefficient Tabulate 


Data 

Point 

No. 

Speed 

rpm 

AP 

in. H 2 0 




Test 

Date 9- 

66 

50 7 

0.28 

21.5 

67 

2,460 

1.29 

22.6 

68 

3,600 

1.84 

23.2 

69 

5,010 

2.46 

24.1 

70 

6,360 

3.04 

24.8 

71 

7,470 

3.39 

26.1 

72 

699 

0.39 

21.2 

73 

1,194 

0.66 

21.4 

74 

1,338 

0.74 

21.4 

75 

1,530 

0.85 

21.3 

76 

1,710 

0.94 

21.5 

77 

1,950 

1.05 

22.0 

78 

2,190 

1.17 

22.1 

79 

270 

0.15 

21.3 

80 

627 

0.35 

21.2 

81 

867 

0.48 

21.4 

82 

762 

0.42 

21.5 

83 

1,374 

0.75 

21.7 

84 

1,647 

0.90 

21. 7 

85 

2,229 

1.19 

22.2 

86 

2,970 

1.54 

22.8 

87 

3,960 

2.00 

23.5 

88 

4,980 

2.45 

24.1 

89 

5,760 

2.78 

24.6 

90 

7,140 

3.35 

25.2 

91 

8,280 

3.80 

25.8 

92 

9,090 

4.12 

26.2 

93 

10,410 

4.63 

26.6 

94 

11,640 

5.12 

26.9 

95 

12,720 

5.57 

27.1 

96 

13,530 

5.95 

27.0 


Viscoaeal No. 1, Gas-Air 


R d 

0 

c in. x 10 


17-69 

9.2 4.18 

44.6 
65.3 

90.9 

115.5 

135.8 

12.7 

21.6 

24.2 

27.7 

31.0 

35.3 

40.0 

4.9 

11.4 

15.7 

13.8 

24.9 

29.8 

40.4 

53.9 

71.8 

90.4 

104.6 

129.8 

150.6 

165.4 

189.5 

212.0 

231.8 

246.7 



- 102 - 


i 


Net Leakage Reduced Data for Seal No. 2, Zero rpm,, Argon Gas 







Data Point No. 

yhg 

P b 

yhg 

1/n k 

N/AP 

,,-14 
x 10 


1 

17,000 

425 

18.00 

17.50 

2 

14,300 

403 

15.10 

15.40 

3 

10,200 

385 

10.90 

11.80 

4 

8,080 

382 

8.68 

9.87 

5 

3,550 

375 

4.02 

6.12 

6 

1,500 

363 

1.91 

4.43 

7 

770 

360 

1.16 

4.05 

8 

2,460 

355 

2.89 

5.26 

9 

5,800 

355 

6.31 

8.13 

10 

440 

3.0 

0.45 

3.81 

11 

145 

1.8 

0.15 

3.68 

12 

251 

2.2 

0.26 

3.81 

13 

680 

4.0 

0.70 

3.92 

14 

65,500 

920 

68.10 

6.76 

15 

40,500 

420 

42.00 

3.97 

16 

28,900 

265 

29.90 

2.93 

17 

64 

1.6 

0.067 

3.58 
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Net Leakage Reduced Data for Seal No. 2, Dynamic Speeds, Argon Gas 


Data 

Point 

No. 

P T 

yhg 

P b 

i/H K 

N/AP 

in -i4 
x 10 

c 

in. x 10^ 



— T 





2500 rpm 

, c x 10 '*= 3.94 in. , r 
* P 

= 212.5 


1 

21,200 

140 

20.70 

13.90 

3.95 

2 

15,700 

101 

15.20 

10.00 

3.91 

3 

11,000 

83 

10.40 

6.88 

3.83 

4 

9,000 

77 

8.62 

4.81 

3.87 

5 

6,200 

70 

5.95 

3.10 

3.87 

6 

4,670 

68 

4.68 

2.35 

4.03 

7 

3,240 

66 

3.27 

1.55 

4.03 

8 

1,850 

65 

1.89 

0.888 

4.03 

9 

29,200 

235 

28.50 

21.10 

3.95 

10 

42,900 

395 

42.40 

33.90 

3.99 

11 

73,000 

978 

72.40 

63.70 

3.99 

12 

1,180 

2.1 

1.11 

1.14 

3.83 

13 

222 

0.6 

0.20 

1.02 

3.71 

14 

97 

0.45 

0.093 

1.017 

3.90 

15 

59 

0.40 

0.055 

1.014 

3.79 

16 

12,450 

53 

12.10 

7.63 

3.95 

17 

2,840 

6.0 

2.67 

1.77 

3.83 

18 

1,410 

2.4 

1.33 

1.15 

3.83 

19 

952 

1.7 

0.92 

1.09 

3.95 

20 

505 

1.0 

0.48 

1.01 

3.83 


30,000 

— 3 

rpm, c x 10 

= 3.51 in. , 

r =246.8 
P 


21 

132,500 

1050 

120.5 

35.7 

3.68 

22 

103,700 

530 

94.1 

16.1 

3.68 

23 

74,000 

150 

62.6 

1.73 

3.44 

24 

88,500 

515 

70.3 

5.19 

3.22 
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Net Leakage Reduced Data for Seal No. 2, Dynamic Speeds, Argon Gas 



Data 

Point 

No. 

P T 

yhg 

P b 

yhg 

1/n k 

N/AP 

in -14 

x 10 

c 

in. x 10 3 


5000 rpm, c x 10 3 = 3.92 in. , r = 891.9 


25 

29,800 

200 

29.40 

15.6 

3.99 

26 

44,500 

360 

43.50 

28.2 

3.95 

27 

78,800 

1000 

77.30 

59.0 

3.95 

28 

4,100 

2.9 

3.61 

0.509 

3.59 

29 

2,850 

1.5 

2.59 

0.352 

3.71 

30 

6 ,660 

8.5 

6.20 

1.25 

3.79 

31 

1,670 

0.75 

1.45 

0.199 

3.54 

32 

525 

0.25 

0.41 

0.111 

3.19 

33 

760 

0.30 

0.65 

0.199 

3.47 

34 

388 

0.20 

0.34 

0.157 

3.59 

35 

355 

0.20 

0.31 

0.163 

3.54 

36 

25,100 

137 

24.70 

11.7 

3.99 

37 

15,800 

81 

- 15.70 

5.96 

4.03 

38 

9,500 

59 

9.26 

2.00 

3.95 

39 

19,100 

106 

18.10 

8.14 

3.97 

40 

12,000 

80 

11.40 

3.64 

3.77 

41 

6,330 

69 

6.30 

0.791 

3.77 

42 

5,200 

75 

4.97 

0.561 

3.75 


10 

,000 rpm, c x 

10 3 = 3.85 in. , 

7 = 1879.2 
P 


43 

17,100 

70 

16.3 

1.11 

3.87 

44 

27,800 

109 

26.8 

5.72 

3.91 

45 

33,500 

135 

31.9 

8.55 

3.86 

46 

23,100 

115 

21.8 

3.23 

3.83 

47 

7,150 

0.6 

6.61 

0.0251 

3.77 

48 

14,600 

5.5 

12.40 

0.331 

3.46 

49 

16,500 

11.0 

14.20 

0.648 

3.51 

50 

10,900 

2.2 

9.71 

0.147 

3.63 

51 

88,100 

1000 

87.10 

52.7 

3.99 

52 

44,400 

265 

42.80 

16.5 

3.91 

53 

36,100 

180 

34.80 

10.7 

3.91 

54 

66,000 

550 

64.50 

33.0 

3.95 



Rarefied Sealing Coefficient Reduced Data for Seal No. 2, Argon Gas 



Data 
Point 
No . 

P T 

Ohg 

P b 

y h g 

1/N k 

A 

R 

e 

c 

c 

in. x 10 3 


5000 rpm, c x 10 3 =» 3.92 in. 


SP1 

4,600 

90 

4.41 

35.0 

0.296 

3.91 

SP2 

35,000 

28,000 

59.80 

22.2 

4.006 

3.94 

SP3 

33,000 

25,800 

55.77 

21.6 

3.739 

3.94 

SPA 

12,000 

4,400 

15.70 

20.0 

1.053 ^ 

3.98 

SP5 

24,000 

16,600 

38.90 

20.6 

2.608 

3.98 

SP6 

8,950 

2,330 

10.80 

23.0 

7.245 

3.98 

SP7 

8,250 

1,670 

9.24 

24.4 

6.195 

3.87 

SP8 

3,180 

50 

2.95 

53.7 

0.198 

3.79 

SP14 

42,000 

34,800 

73.80 

21.0 

4.945 

3.99 

SP15 

56,200 

49,000 

102.80 

20.3 

6.893 

4.06 

SP18 

2,610 

66 

2.36 

70.7 

0.158 

3.66 

SP19 

520 

3.8 

0.45 

364 

0.030 

3.58 


10,000 

rpm, c x 

10 3 = 3. 86 in. 



SP9 

39,000 

23,000 

58.20 

19.7 

7.804 

3.90 

SP10 

14,000 

695 

13.54 

24.6 

1.817 

3.83 

SP11 

10,700 

51 

9.91 

30.8 

1.329 

3.83 

SP12 

19,000 

3,380 

20.80 

20.7 

2.788 

3.86 

SP13 

24,900 

8,160 

31.00 

18.9 

4.162 

3.90 

SP16 

54,000 

36,200 

87.30 

16.8 

11.704 

4.02 

SP22 

4,450 

3 

3.87 

83.1 

0.519 

3.61 


30.00C 

i rpm, c x 

10 3 = 3. 

.23 in. 



SP20 

111,200 

42,200 

123.00 

18.8 

49.463 

3.33 

SP21 

50,200 

1.1 

36.10 

32.1 

14.534 

2.99 

SP17 

63,000 

205 

51.1 

20.3 

20.564 

3.36 



-106- 


Continuum Sealing Coefficient Tabulation Viscoseal No. 2, Gas-Air 



Data 

Point 

Speed 

AP 

A 

R 

c 

o 

No. 

rpm 

in. H 2 0 


c 

in. x 10 J 


1 

543 

0.48 

13.43 

9.73 

4.18 

2 

1,148 

0.99 

13.77 

20.6 

4.18 

3 

4,140 

3.50 

14.31 

73.7 

4.14 

4 

3,480 

2.98 

14.13 

61.9 

4.14 

5 

2,850 

2.48 

13.91 

50.7 

4.14 

6 

2,310 

2.00 

13.98 

41.4 

4.14 

7 

1,710 

1.50 

13.80 

30.4 

4.14 

8 

1,110 

0.99 

13.57 

19.7 

4.14 

9 

570 

0.52 

13.27 

10.11 

4.14 

10 

4,800 

4.01 

14.49 

85.5 

4.14 

11 

5,430 

4.49 

14. 70 

96.6 

4.13 

12 

6,090 

5.00 

14.88 

108.1 

4.12 

13 

6,690 

5.46 

15.04 

118.6 

4.11 

14 

7,380 

6.01 

15.15 

130.6 

4.10 

15 

8,040 

6.51 

15.31 

142.0 

4.09 

16 

8,700 

7.02 

15.51 

153.0 

4.07 

17 

9,300 

7.47 

15.90 

162.0 

4.03 

18 

10,020 

8.01 

16.21 

173.4 

4.00 

19 

10,590 

8.48 

16.51 

181.5 

3.96 

20 

11,250 

9.01 

16.52 

193.0 

3.96 

21 

11,790 

9.46 

16.48 

202.4 

3.96 

22 

12,510 

10.05 

16.46 

214.9 

3.96 

23 

13,110 

10.54 

16.45 

225.3 

3.96 

24 

14,520 

11.74 

15.80 1 

254.3 

4.03 

25 

13,680 

11.07 

15.78 

239.4 

4.03 


Vy 


V 




-107- 


Continuum Sealing 

Coefficient 

Tabulation Viscoseal No. 

2 , Gas-Air 


Data 

Point 

No. 

Speed 

rpm 

A 8 

in.-H 2 0 

A R 

e 

c 

c 3 

in, x 10 J 


Test Date 9-2-70 


26 

4,950 

4.11 

14.30 

89.1 

4.18 

27 

10,200 

8.00 

15.28 

183.6 

4.16 

28 

14,850 

11.83 

15.19 

267.2 

4.14 

29 

16,440 

13.31 

14.95 

296.4 

4.14 

30 

18,180 

15.01 

14.66 

328.5 

4.14 

31 

19,590 

16.54 

14.33 

354.6 

4.14 

32 

21,210 

18.34 

14.00 

384.8 

4.14 

33 

22,650 

20.14 

13.61 

411.9 

4.14 

34 

23,700 

21.67 

13.24 

431.7 

4.14 

35 

24,930 

23.26 

12.97 

455.0 

4.14 

36 

25,920 

24.79 

12.65 

474.0 

4.14 

37 

27,000 

26.55 

12.31 

494.8 

4.14 

38 

27,900 

28.08 

12.02 

512.3 

4.14 

39 

28,890 

29.94 

11.66 

531.7 

4.14 

40 

29,760 

31.81 

11.32 

549.0 

4.14 

41 

30,510 

33.41 

11.05 

563.9 

4.14 

42 

31,410 

35.58 

10.67 

582.1 

4.14 

43 

32,130 

37.32 

10.41 

596.7 

4.14 

44 

32,940 

39.33 

10.12 

. 613.2 

4.14 

45 

34,110 

42.43 

9.72 

637.3 

4.14 



